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Abstract 

Spin chains with exact supersymmetry on finite one-dimensional lattices are considered. 
The supercharges are nilpotent operators on the lattice of dynamical nature: they change 
the number of sites. A local criterion for the nilpotency on periodic lattices is formulated. 
Any of its solutions leads to a supersymmetric spin chain. It is shown that a class of special 
solutions at arbitrary spin gives the lattice equivalents of the N = (2,2) superconformal 
minimal models. The case of spin one is investigated in detail: in particular, it is 
shown that the Fateev-Zamolodchikov chain and its off-critical extension admits a lattice 
supersymmetry for all its coupling constants. Its supersymmetry singlets are thoroughly 
analysed, and a relation between their components and the weighted enumeration of 
alternating sign matrices is conjectured. 
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1 Introduction 

The aim of the present paper is to study models of statistical mechanics, in particular spin 
chains, with exact space-time supersymmetry on the lattice. It is a systematic extension of the 
works on the XXZ and XYZ models investigated in [1] [2] . 

One of our main motivations is that many lattice modelsj^ are known to have continuum 
limits which possess supersymmetry. A rather prominent example is the tricritcal Ising model: 
its scaling limit is described by a A/" = 1 superconformal field theory |5J [5] • A^ = 1 and 

^We focus on two-dimensional lattice models which originate from statistical mechanics, and their one- 
dimensional quantum mechanical pendants (spin chains). Our discussion therefore excludes the wide field of 
lattice gauge theories with supersymmetry, pioneered in [3] (see e.g. |4] for a comprehensive overview). 
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Af — 2 were also observed the Ashkin- Teller model at special values of the coupling constants 
[3 [HI El Hm • Trigonometric vertex models with Uq(sl2)-symmetry were studied in [TTJ [T^l US] • 
The general outcome of these studies is that the supersymmetry occurs at some spin-anisotropy 
commensurable points - an observation which we will encounter in the present work, too. 
Models related to coset theories and their off-critical extensions are studied in and in 

the review paper [T^, and height models in [T7] . 

Even though this variety of examples has been known for a long time, it appears that an 
explicit lattice construction of the supercharges on the lattice has been accomplished in only 
rather few cases. After an early attempt in [TH], the first successful construction were the 
so-called Mg models introduced in [TUJ [20] (an excellent and comprehensive introduction to 
the subject can be found in |21j). They describe spinless fermions on a one-dimensional lattice 
with the exclusion rule that at most £ consecutive sites may be occupied. The models possess 
an explicit supersymmetry whose supercharges insert or take out fermions. In particular, 
the A4i model (fermions with strong repulsion which forbids any two adjacent sites to be 
simultaneously occupied) has been studied for various boundary conditions [22] [23] [23] . Off- 
critical deformations are achieved through staggering of the coupling constants [25] |2H] [27] [25] , 
and their ground states appear to be related to classically integrable hierarchies. Moreover, 
higher-dimensional generalisations of these model were considered in (25] IMl ISU [HI ISS] , and 
display various exciting features such as extensive ground-state entropy and superfrustration, 
and connections to rhombus tilings [31] [33] . 

The topic of the present work is explicit lattice supersymmetry for spin chains. Yet, the 
fermion models will serve as a great source of inspiration. The prime example for a connection 
between the two worlds is the spin— 1/2 XXZ chain with anisotropy A = —1/2 whose continuum 
limit corresponds to a superconformal field theory with central charge c = 1 (a free boson, 
compactified at a special radius). Using a mapping to the Mi model, the lattice supersymmetry 
of the spin chain was discovered by Yang and Fendley [T] . The supercharges have an unusual 
feature: they act non-locally and change the number of sites. Therefore, they are dynamical. It 
is worthwhile stressing that the concept of dynamical symmetries emerged also in the spin-chain 
description of scattering amplitudes and Af = A super Yang-Mills theory [3S] [37] [35]. The XXZ 
results were subsequently extended to the spin— 1/2 XYZ chain along an off-critical extension 
of the A —1/2 point: the lattice supersymmetry was constructed in |2], and shown to be 
manifest within the Bethe-ansatz for the eight- vertex model. It is a quite remarkable feature of 
the spin chains that they possess two sets of supercharges. Their algebraic relations yield a 
lattice counterpart of the Af = (2, 2) supersymmetry algebra. Thus, the findings for the XYZ 
chain constitute a lattice equivalent of the non-local supersymmetry in the sine-Gordon field 
theory at the supersymmetric point |39j. The present article attempts to generalise the concept 
of dynamical lattice supersymmetry to spin chains with higher spin. This requires some general 
tools for the construction of nilpotent operators, and thus supercharges, on the lattice, and 
allows us to find spin chains with lattice Af = (2, 2) supersymmetry at arbitrary spin. 

A very characteristic feature of supersymmetric theories is the special properties of their 
ground states. Indeed, only a short glimpse at the existing literature on both, the XXZ chain 
at A — —1/2 and the A4i model show indeed that their ground-state components display 
various surprising relations with combinatorial problems such as the enumeration of alternating 
sign matrices, plane partitions, fully-packed loops [10] El [H] [22]. In the XXZ-case and 
for the related 0(1) loop model, these observations have pushed the development of rather 
powerful techniques, such as the application of the quantum Knizhnik-Zamolodchikov equation 
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to combinatorial problems |43l Hll H51 H51 147j . and more combinatorial methods using dihedral 
symmetries [3SII12]- This allowed to prove the early conjectures. Off-critical extensions to the 
XYZ model are known, too [301 EB 152] ■ Yet, in many cases a combinatorial interpretation 
is still missing (see however \b'6l I54j ). Despite the omnipresence of lattice supersymmetry in 
these models its relation to combinatorics remains to be understood. The machinery which is 
commonly used to analyse the ground states belongs rather to the world of quantum integrability. 
Yet, once established in a lattice model supersymmetry may be used as a heuristic tool in order 
to identify new interesting (and perhaps even combinatorial) ground states. In this work we 
present the (twisted) Fateev-Zamolodchikov spin chain as an example where this strategy is 
most fruitful. In this case, we determine a lattice equivalent of the Af = (1, 1) supersymmetry, 
which is known to emerge in the scaling limit and does surprisingly not present any obvious 
relation to the fermion models. Guided by this supersymmetry we point out various relations 
of its ground states and the weighted enumeration of alternating sign matrices |55j . 

The layout of this paper is the following. We start with a general considerations in 
section [2] where we review the supersymmetry algebra, and outline the strategy to construct 
representations which are relevant for spin chains with periodic, twisted, and open boundary 
conditions. Section |3] discusses a special class of models for arbitrary spin, the so-called 
trigonometric models which are related to vertex models built from the six-vertex model 
through the fusion procedure. We proceed with a detailed discussion of the spin-one case in 
section [4] and illustrate that the basic requirement for a supercharge to be nilpotent may lead to 
different families of inequivalent model. In particular, we show that the Fateev-Zamolodchikov 
spin chain and its off-critical deformation admit a lattice supersymmetry for arbitrary coupling 
constants. We discuss the zero-energy states of this model, and point out relations to the 
weighted enumeration of alternating sign matrices. Next, we present two other examples of 
supersymmetric spin— 1 chains: on the one hand a spin chain which is closely related to the 
supersymmetric t — J model, for which we obtain an off-critical supersymmetry-preserving 
deformation, and on the other hand the so-called "mod— 3 chain" which interpolates between a 
spin— 1/2 and a spin— 1 model. In section |5] we present our conclusions as well as a variety of 
open problems and directions for further investigations. Some technical details are relegated to 
an appendix. 

Along the main text, we present several results based upon numerical analysis of small 
systems, for example exact diagonalisation. Many of them hint at structures which should be 
present for arbitrary system sizes, and therefore we formulate them as conjectures. 

2 Lattice supersymmetry 

In this section we describe the concept of dynamical supersymmetry for spin chains. We start 
with a brief reminder of the well-known J\f — 2 supersymmetry algebra in section [2?T| In section 
|2.2| we adopt this algebra to the setting of spin chains with periodic and twisted boundary 
conditions, and discuss the general structure of their Hamiltonians as well as basic symmetries. 
Moreover, we present the lattice version of the Af = (2,2) supersymmetry algebra. Open 
boundaries are briefly discussed in section |2.3[ Eventually, we review the example of the XYZ 
chain along its combinatorial line in section |2.4[ 
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2.1 Supersymmetry algebra 



The Af = 2 supersymmetry algebra is generated by two supercharges Q,Q^ , the fermion number 
F and a Hamiltonian H. We denote by H the Hilbert space the algebra acts on. The fermion 
number induces a grading so that 

oo 

H = 0H/, (1) 

where "H/ is the subspace of all states with / fermions. The supercharges are nilpotent mappings 
Q : Hf ~^Hf+i, -.Hf Hf^i which add or remove a fermion. We thus have the relations 

q2 ^ (Qt)2 ^ 0, [F, Q] = Q, [F, Qt] = _gt. 

Finally, the Hamiltonian H : Jif ^ Tif is given as an anticomniutator of the supercharges 

H = Q^Q + QQl 

It commutes with the supercharges Q, and the fermion number F. 

The definition of the Hamiltonian implies that its spectrum is positive definite: any solution 
to H\tp) — E\tp) has E > 0. Each eigenvalue E' > is necessarily doubly degenerate. It is easy 
to see that the corresponding eigenstates organise in pairs such that Q^^/j) — 0. 

They are called superpartners, and their fermion numbers differ by one. Conversely, zero-energy 
states must solve Q\4') = and Q^ip) = Oj therefore do not have superpartners. They are 
called supersymmetry singlets, and are known to be in a one-to-one correspondence with the 
elements of the quotient space ker Q /im Q . 

2.2 Supercharges and spin chains 

Hilbert space. Our aim is to construct representations of the J\f — 2 supersymmetry algebra 
for spin chains. The crucial idea is to identify the number of sites with the fermion number. 
Hence, the supercharges insert or remove a site. At each such site lives a spin-£/2 with 

£ = 1,2,3, We thus have a local quantum space V ~ C^"*"^, and denote its canonical basis 

by |to), TO = 0, . . . , ^ (we use the bra-ket-notation) . The basis vectors are orthonormal with 
respect to the standard scalar product: (TO|n) = 6mn- The Hilbert space of a chain of length N 
is given by the iV-fold tensor product . The canonical choice for its basis consists of taking 
simple tensor products of local basis vectors. Hence, a basis vector is given by a sequence of 
integers toi, TO2, . . . , tojv- We use the common notation 

N 

|toi, TO2, . . . , tun) = |TOj). 

J=l 

The spin is measured by an operator s^: on the basis {|to)} it acts according to s^\m) = 
(m — £/2)\m). Hence, is just the third member of the three su(2) generators at spin-€/2. If 
acting on the j-th factor of V'^^ we write s^ (and do so for any local operator). The total spin 

of any state in V'^^ is measured by Sf^ = ^j- Here, the index N indicates the number 

of sites, what is sometimes necessary as we deal with operators which change the number of 
sites. Sometimes, we will also think of M = 3% + iN/2 as the total ntimber of particles as |to) 
can be interpreted as a state with to particles at a given sites. This will prove to be useful in 
section [31 
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Translation invariance, lattice supercharges. In order to implement the supersymmetry, 
we have to specify the spaces in ([TJ. They are subspaces of V'^^ whose precise form 
depends on the boundary conditions. Most of this work is concerned with periodic or twisted 
boundary conditions (we restrict our considerations to diagonal twists). We denote by (f) the 
twist angle, and introduce a shift operator Tn{(I)) = T/ye"^^" where Tjv is the usual translation 
operator acting on basis states for TV sites according to 

Tjv|mi,m2, . . . ,mAr_i,TOAr) = \mN ,mi,m2, . . . ,mN~i) ■ 

The Hn are related to the eigenspaces of the twisted shift operators. We will sometimes call 
them (pseudo)momentum sectors. In order to single out the relevant sectors, we need to explain 
the construction of the supercharges. These are built from local operations : V^-^ — >■ y^C^+i) 
which are related by translations. They are defined as 



(-1) 



V 



q (g) 1 (g) • • • «) 1 , J = 1, . . . , iV, 
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and therefore transform a spin at site j to a pair at sites j,j + 1. The matrix elements for this 
transformation are encoded in the local supercharge q : V ^ V V. Moreover, the operation 
is weighted by a site-dependent string (—1)-'^^ which will be crucial to make the supercharge 
nilpotent. We find 

qj+i = -TN+iiq^)q^TN{c^)-\ j = 1, . . . , iV - 1. 

On a system with N + 1 sites, a pair has + 1 possible positions. This motivates the 
introduction of another local operator qo through extension of the preceding equation to j = 

qo = -W(0)-iqiT^(0) = (-l)^T^^+l(0)q^^. 

This operator acts always on the last site (irrespectively of the number of sites) . We see that 
there are two definitions for qg: (i) shifting the action of qi to the left or (ii) shifting the action 
of qAT to the right. These definitions have to be compatible what is true only if the sign (—1)^ 
is present in the second definition. Because of this requirement we find that the 'H n are the 
eigenspaces of Tm{4)) in V®^ with eigenvalues — (—1)^+^. For a non-vanishing twist angle, 
this imposes a constraint on the magnetisation. Indeed, it follows from TN{(j))^ — e"^'^" that 
we need to restrict our considerations to subsectors where 

4>S% = mod 2tt. (2) 

Often, a non-zero twist angle cannot be chosen arbitrarily. There are constraints coming from 
the local operator q: its components, defined through 

e 

have to solve the equation 

a„.,,k (e'*(^/2+"-^-^-) - l) = 0. (3) 
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Once again, this is a consequence of the two definitions of qo given above. This equation holds 
always for (f) — 0. Below, we will see that non-zero (/)'s are often quantised: only very particular 
values are compatible with this condition. 

With all these preliminary definitions we are ready to define the global supercharge acting 
on 'Hat as the following sum of the local operators: 




In all other (pseudo)momentum sectors we set Qn = 0. It is then straightforward to see that 
Tpfj^i{(f))QfqTjq{(f))^^ — —Qn on "H^r. We thus see that Q^q does not only insert a site but also 
increases the (pseudo)momentum by tt. It is worthwhile stressing that the construction of the 
supercharges, the restriction to special quantum sectors as well as the quantisation of twist 
angles fits well into the general framework of non-local currents in lattice models as developed 
by Bernard and Felder '57] . 



Nilpotent operators. Our aim is to choose the local operator q -.V^V^Vso that Qm 
is nilpotent in the sense that 

Qn+iQn = 0. (4) 



We give two very simple criteria for this property to hold: (i) the proposition 2.1 rules out 
redundant parameters in q and leads to the concept of gauge-equivalent local supercharges; ( ii) 
a criterion for the global supercharge to be nilpotent is given in proposition |2.2[ In both cases, 
the arguments are based on simple telescopic cancellations. 

Proposition 2.1 Suppose that there exists a\x) & V such that q|m) = |m) ® |x) + |x) ® \m) 
for every \m) g V . Then the global supercharge Qn vanishes. 

Proof: It is sufficient to verify the statement on basis vectors. Let us first consider a simple 
spin configuration \fi) = \mi, . . . ,mi\[) for any admissible choice of the labels rrij. Then 

-I- |toi, . . . ,mi_i) (g) Ix) (8) \mi,...,mN)) 

and 

qolA^) - - {(TNi^M) ® \x) + \X) ® Im» • 
The sum over all q^ leads to telescopic cancellations, and we are left with 

Eq, I Im) = ((-i)^+V> ® lx>- 

From any state we can construct an eigenvector of the shift operator, given by 

N-l 
J=0 
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Using the preceding formula, we find that 
which vanishes for t^q = (—1)^+-'^. □ 

We conclude that any q contains £ + 1 redundant parameters as long as no other restrictions 
are imposed. Moreover, we see that two local supercharges q, q' are gauge-equivalent in the 
sense that if 

q' = q + qx' qxl™> = l*^) ® lx> + lx> ® |™), 

then the global supercharges coincide Qjsi — Q'j^. This provides a way to remove redundancies. 
Once this is done we need a more refined criterion in order to decide if Q holds. Crucial for 
the property of Qn to be nilpotent is the implicit string (—1)^^^ in the definition of q^. It 
leads to a set of anticommutation relations: 

q^qj + qj+iqi = o, i < ^ < j < a^, (Sa) 
qoqj + qj+iqo = 0, j i,...,iv- i. (5b) 

They simplify considerably the expression Qn+iQn, and leave only a few terms to be analysed. 
Guided by the nature of the preceding argument, which is based on telescopic cancellations, we 
are led to the following 

Proposition 2.2 The supercharge is nilpotent if there exists a vector |x) G such that 
((q (g) l)q - (1 (g) q)q)|m) = Ix) «) - Ix), for all \m) e V. (6) 



Proof: Notice that because of the anticommutation relations (|5| we can reduce the expression 
to 




We investigate the action of the different terms on a simple basis vector |/Lt) = |mi, . . . , mjv). 
For i = 1, . . . ,N we have 

(qi+i + qi)qi|Ai> = l^i, . . . ,mi)(g)|x) <E) \m,+i, . . . ,tojv) 

- |mi, . . . ,mi_i) (g) Ix) (g • ■ -jruN). 

The remaining two terms are slightly different: 

(qi +qo)qo|A*> = ((T^(0)|m» -g |x» - Im> ® lx> 

(q^+iqo + qoqjv)|A^) = (-1)"+^ ((Tiv(^)|M» ® Ix) " Tjv+2(</>) <8> |x») • 
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Hence, we find 



J2 (qj+i + qj) qj +qjv+iqo + qoqA^ I Im) =Tjv+2(0) {{TNm^i) - ® \x)) 

+ ((-i)^+it^(0)|m)-|a^))®Ix). 

From we can construct a momentum state like in the preceding proposition. As Qn 
only acts on states with TAr(0)|V') = {—l-)^~^^\fp) it follows that Qn+iQn\''P) — O 

The equation ^ provides an eSicient algorithm to find supercharges at spin— £/2. The 
general strategy is the following. First, choose a rectangular + 1)^ x {i + 1) matrix q 
with possible restrictions (some entries might be zero due to desired symmetries such as 
spin-reversal symmetry, particle number conservation etc.). Next, determine \x) by solving the 
linear equations for its components, defined through \x) = n=o Xmn\'m, n). In terms of the 
components of the local supercharge am,ij we find two expressions 

/ \ / I. ^ ) 

Z^j—Q \^k,jn^j^kni ^k^kj^j^mn) 7 Ti^ K. 

These expressions must hold for any allowed value of fc. Third, one is generically left with a 
number of quadratic equations for the non-zero components of q whose solutions lead to a 
nilpotent supercharge. There may be many solutions, some of them inequivalent in the sense 
that q and q' cannot be related through q' = (u (g) M)qu~^ (where u is a unitary transformation 
on V) and/or a gauge transformation. It would certainly be interesting to classify them. In 
the following chapters, we give various interesting examples. 

Hamiltonian. Given a solution to (|4]), we construct the Hamiltonian as the anticommutator 

-H^w = Q]vQw + Qw-iQIv-i- (^) 

Notice that its action is non-zero only on the subsectors where the translation operator T/v(0) 
has eigenvalue (—1)^"''^. On these sectors, we have all the features of a supersymmetric 
model explained in section |2.1[ The spectrum of the Hamiltonian is positive. Possible zero- 



energy eigenstates (supersymmetry singlets) are annihilated by both Qjv and Qjy-i' 
in one-to-one correspondence with the quotient space = ker Qn /imQff_i. All other 
eigenstates are organised in supersymmetry doublets dV'), QatIV'))- due to the property that 
Hn+iQn = QnHj^ the superpartners have the same energy. 

For later developments it will be useful to understand how the Hamiltonian Hn acts locally. 
In fact, the equation (|8| can be simplified as there exists a second set of anticommutation 
relations, this time between the local q^ and their adjoints: 

q,q] + qj+iq, = 0, 1 < i < J - 1 < iV - 1, 
qoq]+q]+iqo = 0, 2 < i < N ^ 1. 
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As before, these relations are a simple consequence of the shift and the implicit string in the 
definition of the q^, and do not depend on the actual structure of q. Using these, we find that 
the Hamiltonian can be written as sum over nearest-neighbour interactions: 




/i,- ,,+1 n 



TV- 



Here, Hn — N^^J2j=o l)''^^^'-'7Ar(0)'' is the projector on the momentum spaces we 
are interested in, and ft-jj+i is the Hamiltonian density acting on sites j,j + 1. We have 
/ij+ij_l_2 = T]\[((f>)hj,j^iTff((f>)~^ with the site identification + 1 = 1. As a matrix, the 
density is given by 

h = -(q^ (g) 1)(1 (g q) - (1 (g) q''')(q g) 1) + qq'*' + ^ (q'^'q g) 1 + 1 g) qU) ■ (9) 

This choice is of course not unique because of translation invariance, yet it appears natural 
because h becomes self-adjoint. 

Given that the local supercharges q are not unique we ask how the Hamiltonian density 
h defined in (|9| changes if we deform q. Proposition 2.1 provides a criterion for redundant 



parameters within the local supercharges: if we have q' = q + q^ with qxI'Ti) = |x) (g 
I™) ® \x), \x) G then Q^r = Q'n- Yet the densities h and h' differ in general. Indeed we 
find 

/i' = /i + (A (g 1 - 1 (g A) , A = ^^Xj{ai,jk - tti^kj + akji - ak,ij)\i){j\- (10) 

ijk 

The term A g) 1 — 1 g) A deserves the name 'local botmdary' for obvious reasons. We see that 
it vanishes obviously if atjk = cii^kj- When summing over all sites the contributions from the 
local boundary terms cancel out telescopically and the total Hamiltonians for q and q' are the 
same. Yet, local boundary terms may be helpful in order to understand local symmetries. We 
will see this when studying quantum group symmetries in sections |3.3| and |4.1.3| 



Basic symmetries. If the local supercharge q has some symmetries, then they are naturally 
inherited by the Hamiltonian density h, and may lead to symmetries of the Hamiltonian H^. 
We will be particularly interested in cases related to particle number conservation/conservation 
of the magnetisation. We thus analyse if there are suitable a G (0, 2tt) such that the equation 

holds. We distinguish two special cases. A supercharge is called trigonometric if this equation 
holds for arbitrary value of a in case of which the total magnetisation is conserved. This is the 
case if (s^ g) 1 + 1 g) s^)q = q(s^ -I- m) for a fixed number m S Z/2. It is called elliptic if the 
condition on h holds for a = tt. In that case the total magnetisation is conserved mod 2. The 
denomination is clearly inspired from the context of integrable vertex models whose i?-matrices 
have generically trigonometric weights if some particle number conservation is imposed, whereas 
the weights are expressed in terms of elliptic functions if this conservation holds only mod 2. 
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However, these are not the only interesting cases. Indeed, in section [43] we discuss a case where 
a = 27r/3 so that the magnetisation is conserved mod 3. 



We discuss briefly two other symmetries for closed chains with periodic boundary conditions 
and open chains. We start with the spin-reversal operator defined through 

i?|m) ^\e-m). 

If q is invariant under spin reversal in the sense that (R (g) R)qR = ±q then h will be invariant 
[R (g, R^h] = 0, and thus [Rn,Hn] = where Rn = i?^^. While the local mvariance 
under spin reversal remains valid in the case of non-zero twist angle, the global symmetry 
operator needs to be replaced by R'j^ = RnK where K denotes complex conjugation. The 
simplicity of all these statements is due to the fact that we deal with a local operation. 
This is not the case for the reflection operator acting on a state for N spins according to 
PN\mi,m2, ■ ■ ■ ,mjv-i,TOjv) = \mN,mN-i, ■ ■ ■ ,TO2,"t-i)- We find 

PN+iqj - (-l)^+iq^%+iPjv, J = 1, . . . ,iV, F^+iqo - q^PNTN- 

The local supercharge q°P differs from the action of q by a transposition: q°P|m) = J2jk'^m.jk\kj). 
The global supercharge Q n will have a definite transformation behaviour with respect to reflec- 
tion symmetry if q°P = eq with e = ±1 what is equivalent to am,jk = ^cim,kj- Then it follows 
that Pnj^iQmPn = (— l)^^^e Qm on T-Lm what implies the invariance of the Hamiltonian under 
reflection [P^Vi^^Af] — 0. 

M = (2, 2) lattice supersymmetry. In some cases it is possible extend the preceding 
considerations by introducing a second set of supercharges Qat, QJv- They are subject to the 
following anticommutation relations: 

QnQn-1 = QnQn-1 — 0, QIv-iQIv — Qjv-iQ]v 

QnQn + Qn-iQn^i — 0, QnQn + Qjv-iQIv-i 



QnQn-1 + QnQn-1 — 0, QIv-iQat + Q^N-lQ^ 



N 



0, 


(11a) 


0, 


(lib) 


0. 


(He) 



To each pair of supercharges is associated a Hamiltonian 



Hn — QnQn + Qn-iQ^n^i, — Q^j^Qn + Q^^iQ^j^^^. (Hd) 

This set of relations is a lattice version of the J\f — (2, 2) supersymmetry algebra known from 
field theory (see e.g. [35]), here with zero centre. In the field theory, non-zero central extensions 
occur on non-compact spaces [55] ■ We do not expect them in the present setting of finite 
lattices. Moreover, in the field theory the operators Q and Q are usually associated with left- 
and right-movers, a concept which is absent on the lattice. Here, we mean by A/" = (2,2) merely 
the presence of a second copy of the supersymmetry algebra. 

From the relations ( [TT] ), it follows that (resp. Hn) commutes with Qn, Qjy (resp. 
Qn, Q\\f), and furthermore [Hn , Hn]=0. Thus, they can be diagonalised simultaneously. The 
eigenvectors to strictly positive eigenvalues organise in quadruplets with one state at — 1 
sites, two at N sites, and one at -I- 1 sites: 

m,QN-i\^),QN-iW,QNQN-im, Qjv-alV') = Q^N-2W = 0- 
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The two states QAr_i|?/'), Qn -il^P) in the middle of the quadruplet can be mapped onto each 
other by the nilpotent operator 

Cn = QnQn, Cff = 0. 



It commutes which both H^, H^q as a consequence of (111 and is generically a rather non- trivial 
(non-local) conserved charge. We will be interested in the case where the two supercharges are 
related by a similarity transformation Un acting on Hn: Qn — Un+iQnUJ^^ . This implies 
that the two Hamiltonians are conjugate Hn — UnHj^U^^ . In section [s] we will encounter 
a case, where Um is a symmetry of the Hamiltonian, so that _ffjv = but does not leave 
invariant the supercharges. 



2.3 Open boundary conditions 

Consider a q such that ^ holds with |x) 0: (q (g) 1 — 1 (g) q)q = 0. This means that the 
mechanism which makes Qm nilpotent is completely local. In this case there is an extension to 
open boundary conditions which has a fully supersymmetric Hamiltonian (no restriction to 
special subsectors in necessary). In the case of an open interval with sites j = 1,2, ... ,N we 
define the global supercharge as 

N 

which is nilpotent if |x) — 0. As before the Hamiltonian is given by an anticommutator: 

N-l ^ 

Hn = QnQn + Qn-iQn-i = 51 '^JJ+i + 2 ("^i*^^ ^ qjyqjv) • 

where /ijj+i is the Hamiltonian density defined in ([9|. The boundary terms can be interpreted 
as surface magnetic fields (this is for example the case for the spin- 1/2 XXZ chain [T]). 



2.4 An example: the XYZ chain 

We review the example of the XYZ spin chain at the combinatorial line [33 ISHl EI], 
discussed in [5] from the point of view of supersymmetry, within the present formalism. The 
supersymmetry is a generalisation of the XXZ-case which acts locally like — )■ and 1 — >■ 00 
pp. The XYZ chain is an elliptic chain in the sense that it breaks particle number conservation 
mod 2. As the XXZ supercharge changes magnetisation by AS^ ~ —3/2 we thus seek for a 
supercharge which allows AS^ = —3/2, 1/2. Thus we are led to a q defined by 

q|0> = aoaollO) + ao,oi|01), q|l> = ai,oo|00) + ai,ii|ll)- 

It turns out that the only interesting choice leading to a nilpotent supercharge is ao,oi = ao,iOi 
i.e. q = q°P. Scaling allows to set Ci^qo — 1. Furthermore, we parametrise ao,oi = A and 



fli ii = — + 2A. Applying proposition 2.1, we see that the A-dependent part does not 
contribute to the global supercharge, and thus we can set A = 0. Next, one checks that ([6| 
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holds with Ix) = — C|00)- Thus, we find bona fide supercharge on the lattice. It generates the 
XYZ-Hamiltonian 



- 4 E (1 + + (1 - o-M,, + + -^^^^ (12) 

where tr", a = 1,2,3 denote the usual Pauli matrices. The Hamiltonian is invariant under spin 
reversal, and thus second local supercharge q — {R ^ R)qR may be defined. Thus, we have 

q|0)=0, q|l) = |00)-C|ll), q|0) = |11)-C|00), q|l) = 0. 

It was shown in [2\ that the corresponding global supercharges generate the Af = (2,2) 
supersymnietry algebra on the lattice. All this holds only for periodic boundary conditions as 
for non-zero ^ the condition (|3| allows only the twist angle = 0. 



Rotational invariance and duality transformation. We use the occasion to inspect a 
known duality transformation of the Hamiltonian (12 1 in parameter space from the point of 
view of the supersymmetry algebra. We denote by p°'{9) — exp {iOa"" /2) a local rotation about 
the a— axis by an angle 9. A global rotation on 7i]y then corresponds to ^%{d) = Y[j=i P"j{^)- 
Moreover, we write q — q(C) in order to stress the explicit dependence on Q. Let us start with 
a = 3. We find that 

{p^e) ® p^e)) q(C)p'(-e) = e^"\{e-^''C:). 

We see that the transformation C — > —C, corresponds to a global rotation about the z— axis with 
angle 6 = 7r/2. Next, let us consider rotations about the j/— axis. Choosing the same angle we 
find 

(p^ W2) « P^(V2)) q(C)P^(-/2) = (^) '-^^^^ - (^) q.e. 



with 



C + 3 



Here qrcd is redundant in the sense of proposition 2.1 with vector \x) = (|0) + |l))/\/2. Therefore, 



it does not contribute to the global supercharge and therefore can be discarded. Thus we have 

nl,,[±./2)Q,[Onl[±./2)-^ = (i^) Q^(C')^Qa^(0 _ 

Hence these rotations imply simple transformations for the supercharges, what is of course 
compatible with the known invariance properties of the Hamiltonian [61) . The strategy to find 
invariances is generic and can be applied to other cases. We present an example for a spin— 1 
chain in section |4l 
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Zero-energy states. The supersymmetry singlets of the XYZ Hamiltonian ( 12 1 were studied 
in |Sni EIJ [23 [52]. These investigations rely on the following existence conjecture: if the 
number of sites A'^ = 2n + 1 is odd then there are exactly two linearly independent zero-energy 
states, whereas no such states exist if A'' = 2n is even. The statement holds for periodic 
boundary conditions. Here we give a short proof of this conjecture in the momentum sectors 
where the supersymmetry exists. It is based on Witten's conjugation argument, which states 
the following. Consider an invertible transformation M onH and let Q = MQM~^ where Q 
is a supercharge, i.e. — 0. Then the Hamiltonians H = {Q,Q^} and H = {Q,Q^} have 
the same number of zero-energy eigenstates. If M leaves invariant certain quantum numbers, 
such as the fermion number, then the result holds sector- wise. We refer to for the details 
of the proof. The basic idea is to exploit that the zero-energy states of H are in one-to-one 
correspondence with the elements of the quotient space = kerQ/imQ, in other words the 
number of linearly independent ground states is equal to dimi^g. Using the conjugation by 
M one establishes a bijcction between the elements of i^g and S)q. This leads to the equality 
of their dimensions, and thus to the desired result. Notice that this does not imply that the 
singlets of H can be obtained from the singlets of H through the transformation M however. 

We apply this argument to the XYZ spin chain by constructing a simple transformation 
M for its supercharge. To this end, notice that for any real number n the two local 
supercharges q(C) and q(^^C) can be related through conjugation as follows 

(m(Ai)®m(/i))q(C)m(M)-^ = M"^q(M'C), with m(M) = |0) (0| + /i|l) (1|. 
Thus, it follows that with A/jv(/^) — ■m{ii)'^'^ we have the relation 

QAr(c,M) - MN+l{^i)QN{C)MNi^^)-' = ^i-'QN{^^\)■ 

These new supercharges generate the Hamiltonian 7Jjv(C:M) = f^~'^Hff{fi'^(). According to 
Witten's argument, it has the same number of zero-energy states as H]s[{(). This statement 
holds for any 0. We suppose now that C > without loss of generahty. Set /x = l/v^i and 
observe that 

N 

hn{c, i/Tc) = (H^ii) = cE(i - 

is, up to a multiplicative factor, the Hamiltonian for a one-dimensional classical Ising model. 
We introduce the states |-l>) = (|0) + |1))/V2, |^) = (|0) - |l))/\/2. In terms of these, the 
only zero-energy states are the fully-polarised configurations |-;— • • • ) and |— >■— >■ • • • — >■). 
Restricting to our momentum sectors with = (—1)^^^, wc find that Hn{1) has two zero- 
energy states for A^ odd, and none for A^ even. The same holds thus for Hn{() for any C > 0. 
Using the duality transformations, it is easy to see that this result extends to any C. 



3 Trigonometric models 

In this section, we present a class of spin— ^/2 chains with dynamical lattice supersymmetry. 
The construction is inspired from the A^^-models for lattice fermions [501 [H]. They possess a 
(non-dynamic) Af = 2 supersymmetry on the lattice. The crucial observation here is that once 
reformulated in the language of spin chains one may impose spin-reversal symmetry, which 
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has no obvious counterpart on the fcrmion side, and thus enforce the existence of a A/" = (2, 2) 
lattice supersymmetry in the sense of (111. In this way we arrive at lattice models for the 
Af = (2, 2) superconformal minimal series with explicit supercharges on the lattice [IH]. 

We proceed as follows. First, we review briefly the definition of the models in section 
|3.1| and discuss the supercharges for the related spin chains. In section [X2] we analyse the effect 
of spin-reversal symmetry, and show that all coupling constants of the theory are essentially 
fixed by this requirement. Next, we prove that suitably modified versions of the corresponding 
Hamiltonians possess a local quantum group symmetry. This leads us to an identification of 
the models with spin— £/2 chains constructed from the six- vertex model through the fusion 
procedure, with a particular anisotropy depending on £. In the last section, we address the 
computation of the Witten index, which indicates the existence of zero-energy states for the 
spin chains if periodic boundary conditions are imposed. 



3.1 Supercharges 

The A4i describe spinless fermions with the constraint that connected fermion clusters cannot 
contain more than £ particles. Their supercharges take out single fermions from a state. This 
operation is weighted by an amplitude am,k: k — 0, . . . ,m — 1 if the (fc + l)-th member in a 
string of m consecutive particles is removed. See figure [ija) for an illustration. There is a 
natural identification of a string of m fermions between two empty sites with a basis vector |m) 
in the vector space V ~ C^+^ as illustrated in figure lljb). We translate the splitting process 



(a) 



-fc-1 



(b) 5 — • — • • — c < — > \m) 

V ^ ' 

m fermions 

Figure 1: (a) The local action of the supercharge in the fermion model: an occupied site is 
represented as •, whereas an empty site corresponds to o. The supercharge splits up a string 
of m particles into two strings with k and m — k — 1 particles, (b) Correspondence between 
string of m fermions and spin state |m). 

shown on figure [ija) into the spin language: the local supercharge takes thus the basis vector 
|m) £ to a pair |fc, to — fc — 1) e V with amplitude am,k- We set Qm.k = for fc < and 
fc > m — 1. Thus we have a local q with amplitudes am.jk — Sj+k+i,mO,m,j'- 

ni— 1 

cilm) ^ ^ ara,k\k,m - k - 1). (13) 

k=0 
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Notice that q changes the magnetisation by + 2)/2. Combining this with ^ we conclude 
that the set of admissible twist angles is given by 

, 47rm f^+l, ^odd, 

1 + 1 1*^/2, £ even. 

Next, we want the global supercharge Qjv built from q to be nilpotent. Hence we need to find 
a vector |x) such that ([6]) holds. Writing out the components in ([t]) with am^jk = Sj+k+i,m<im,j 
we find |x) = 0. It is thus sufficient to solve the equation (1 (g) q)q = (q l)q. In order to 
simplify the discussion, we impose another condition am,k — am,m-i-fe- Hence the splitting 
depends only on the length of the subsequences but not on their order. This requirement 
implies q = q°P and therefore the Hamiltonian will automatically be parity symmetric. Using 
this we obtain the following recurrence relation: 

We set fc = and find 

)rn~n—l / x 
am-l,n =11 an+1,0- (15) 

Hence, it is sufficient to know the numbers Ofe^o, fc = 1, . . . , This is essentially the same result 
as for the M.i models |19]: up to a global factor there are t—\ free parameters. 

3.2 Hamiltonian and spin-reversal symmetry 

We fix the free parameters by imposing spin-reversal symmetry m ^ i — m on the Hamiltonian 
density h. This has no obvious equivalent in the fermion models as the exclusion rule does 
not admit a particle- hole symmetry (at least not for fixed length) . For the matrix elements 

h{r,s),(m.n) = {r, s\h\m,n) this implies 

^(r,s),(m,n) = h(^e^r,i~s) ,{£~mJ-n) ■ (16) 

It is straightforward to deduce the matrix elements in terms of the amplitudes am^k- 

\k=a fe=0 




Our aim is to find the numbers a„j „ such that (16 1 holds. The structure of the Hamiltonian 
density implies that we need to solve a set of quadratic recursion relations what is done in 
detail in appendix We express the result in terms of g-integers defined through 



q-q 1 



We find that supercharge defined through (13 1 leads to a spin-reversal invariant Hamiltonian 
density if we choose the coefficients 



[m — n] [n -|- 1] ' 



m = 1, . . . , n = 1, . . . , m — 1, (17) 
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where the g-integers are evaluated a.t q — exp(i7r/(£+ 2)). Let us sketch how this statement 
is verified for the matrix elements of h. We start with the diagonal 6„i,n — h(^m.n),{m,n)- We 
consider the difference 

I _7 2 _2 |2 |2_2 

Ol-m,l-n "ni," — a2£-(m+ji) + l,£-m '^m+n+l,m + '^m+1,0 "i" °'1,0- 

The construction given in appendix [X] implies that expression is zero for m -\- n = I. We focus 
on m + n < £ (the case m + n > £ can be obtained through symmetry) : the first term on the 
right-hand side vanishes because 2^ — {m + n) + 1 > t + 1. Hence we need to check only the 
equation a^n+i.o + "^^n+i.o ^ '^i o = '^m+n+i m- Indeed, it holds what can be shown from simple 
properties of the integers. Thus, be-rn,i-n = ^m,n- Next, let we turn to the off-diagonal matrix 
elements for which (r, s) ^ (m, n). Without loss of generality, we suppose that r < m. Using 
the explicit form of the matrix elements /i(r,s).(m,n) given above, the spin- reversal symmetry is 
equivalent to 

— 0'2l-{m+n) + l,l-rO-2t-(m+n) + ld-m ~ <^i-r,t-mO'l-n,l+r-{m+n) ■ 

It is obvious that relation holds ii m + n = £. For symmetry reasons, it is sufficient to consider 
thus m + n < £. In this case the first term of the second line vanishes. Moreover, using 
[ni + l]ai-r/-m — [r+ l]am,r, and some basic identities for g-integers the equality follows. This 



concludes the proof of spin-reversal invariance for the constants (17 1. 

In order to make it more explicit, we expand the Hamiltonian density h obtained from the 
trigonometric supercharge in components as follows: 

£ min(m2,£— mi) 

h= ^ ^ I3^^^^jmi + n,m2-n){mi,m2\. (18a) 

mi,m-2=0 n=— inin(mi 7712) 

For n > the coefficients are given by 



Pm„m, j^j y j^^^ ^ ^jj^^^ + n + 1] ' ^ 

where Mi = min(mi,£ — 7712) and M2 — min(m2,i! — mi). The case n < can easily be 
obtained from the symmetry of h: 

If n = we have 

1 ^ [n -f 1] — [n — 1] /ioj\ 

fc=l ^ ' 

The constants Cm appear to be some q-analogues of the harmonic numbers, enjoying the property 
Cm = Q+i-m- Notice the similarity of these matrix elements with those of the spin— i?/2 XXX 
Hamiltonian studied in 621. 
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TV = (2,2) supersymmetry. Imposing the symmetry under spin reversal leads in a natm'al 
way to the question if a second copy of the supersymmetry algebra, generated from q = 
{R (S) R)qR, leads to the Af = (2, 2) algebra. This is certainly the case for £ = 1. Indeed, for 
any ^ > 1 we find indeed that the choice (17 1 for the coefficients implies (11 1. It is sufficient 
to verify the relations locally. This is very similar to proposition |2.2[ first one checks that 
(q (g) 1 — 1 (8) q)q + (q 1 — 1 ® q)q = for the coefficients ( 17 1 . Second, one applies the argument 
of telescopic cancellations in order to prove the relations on the momentum spaces of interest. 

Examples. We illustrate the cases £ = 1 and £ ^ 2. The local Hamiltonian is obtained from 



• For ^ : 



1 we find 
/I 



V 



4' 



1 / 



where we indicated only non-zero matrix elements. This is just the local Hamiltonian for 
the spin- 1/2 XXZ chain at A = -1/2. 

• For ^ = 2, we obtain the 9x9 matrix 
/ 2 



h = 



V2 



\ 



1 


-V2 


-1 


-1 


2 










2 








2 


-1 


-1 


-V2 


1 



-1 



2 / 



This expression reminds strongly of the local Hamiltonian for the Fateev-Zamolodchikov 
chain [53]. The latter is given by 



hpZ - ^ (s" ® S"^ + (S°)' ® 1 + 1 ® (S")2) - J2 ^abS^'s'' ® s'' 



(19) 



a—1 a,fc— 1 

with coupling constants Jq and a symmetric matrix Aat — Aba such that Aaa — J a, and 

Ji = J2 = 1, J3 = cos 277 and A12 = 1, v4i3 = A23 = 2 cos 77 - 1. (20) 

The spin operators s°, a — 1,2,3 are the usual spin— 1 representation of su(2). We 
identify 77 = 7r/4 as a good candidate, but the matrices h and /ipz do not quite coincide. 
In fact, the difference is just a simple local gauge transformation: 



V2 



(21) 
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Therefore, we conclude that for spin one the supersymmetric Haniiltonian is the twisted 
Fateev-Zamolodchikov Hamiltonian (in the sense that boundary conditions are s^_^_^ = 
(— 1)^5]*^, s^_i_j^ = sfjthus antiperiodic for N odd) at 77 = 7r/4. We shall see later in 
section |4] that this is not the only supersymmetry in the Fateev-Zamolodchikov chain. 

3.3 Local quantum group symmetry 

We proceed with a more detailed analysis of the Hamiltonian density: we show that a suitably 
modified version of h commutes with the quantum group U|j(sl2) with q = exp (i7r/(^ + 2)). For 
generic q = exp(i77) the algebraic relations of its generators are given by 

sm rj 

A highest-weight spin— £/2 representation acting on a vector space V ~ C^+^ is given by the 
following action on the basis vectors |rn), m = 0, . . . 

5^\m) = {m-e/2)\m), 

S^\m) — — TO][m -t- l]|m -I- 1), s~\m) — — m + l][m]\m — 1). 

In order to obtain the action of Ug(sl2) on two sites we use its Hopf-algebra structure (see e.g. 
[01]). The comultiplication is defined on the generators according to 

A(s±) =5=^ «)S=^, A{s^) = (g) 1 + 1(g) . (22) 

These few definitions are sufficient to prove the hidden quantum group invariance of h. The 
relations become more transparent if we consider a slightly modified (twisted) version h'. To 
define it , we need the unitary transformation u — 1 g) exp Ittb^ . We set 

h' = uhu^^ + A(s-^ (g) 1 - 1 (g)S^), 

where A is a constant to be adjusted. We see that h' decomposes into a local bulk term uhu~^, 
and a local boundary term A(s'^ (g) 1 — 1 (g) s^). The density h' commutes with Ug(sl2) at 
q = exp {m/{£ + 2)), i.e. [h' , A(s=^)] = [h', A{s^)] = 0, provided that we choose 



A = isin^— j. (23) 

The proof is elementary, and we only sketch it. The commutation relation with A(s'^) is trivial. 
We focus on [h' , A(s+)] = 0. This amounts to prove the identity 

7^+i,„9'/'- V[™ + + 2] + ^itn+ir-'^Win + l][n + 2] 



= iLnl'^'^'-^Virri + J + l][m + J + 2] + litnQ'^+'^'-'^Win - j] [n-J + l], 

m2d-mi) n I 

- inin(mi,£— m2) 'mi, "12 I 



where the are defined through h' = EL^.^^^o -7„'„VT,'^-m2) T'mi,™^ + ™2 



n) (mi , 7712 1 and given in terms of the components of h by 

linn = (-iF^mn + ^5]fl{m - Tl) . 
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We start with j — and fix the value of A. For m + n < £ we find 



Using the simple identity 

[m + n] =q™[n] + g-"[m], (24) 
the right-hand side simplifies to g'"+^/[m + 2]. The left-hand side is given by A -I- ([m -f 3] — 



[m + 1]) /{2[m + 2]). Equating both sites leads to (23 1. The case m + n> £ is treated similarly, 
and gives of course the same result. The constant A appears also if we set j = —1. Yet this case 
is equal to the one just considered because of symmetry of the coefficients under j — >■ — j 
(which they inherit from the numbers Pmn)- Next, we have to verify the equation for j > 1. 
In this case, we have 7^„„ = (— 1)''/3^„. Grouping together the terms containing on the one 
hand l3^^i „, /3^„ and Pt^n+i^ Prim the other, the identity is readily verified by using the 



explicit expressions (18 1 as well as the identity (24 1. As before, the work for j < — 1 can be 
reduced to this case by using the symmetry of the coefficients under j — > — j. Finally, the proof 
of [h' , A(s~)] = follows the same lines. □ 

3.4 Relation to higher-spin XXZ chains 



Looking at the two examples £ = 1 and £ = 2 given at the end of section 3.2 it seems natural 
that the generic case leads to spin— ^/2 versions of the XXZ spin chain. The Hamiltonian 
densities of these spin chains with arbitrary anisotropy were computed explicitly in |65l . Let 
C — s^s+ + [s^][s^ + 1] be the quadratic Casimir of Uqish). Define — A(C), acting on 



y (g) y ~ (g) V.^^^ . In terms of the generators it is given by 

^ «) + (K) S+ H ^ f cos 1] f q"^"" ® q"" + g"" ® q"''''') 

4 sin 77 V V / 

-cos(^+ 1)7; {q^'' (g) g'' (g) 9"''')) 



The Hamiltonian density of the higher-spin XXZ chain is a polynomial of order £ in X^. 
Explicitly it is given by 



'XXZ« 



2E^. n 



2Xi - H[m+ 1] 



ml [7 + TO -|- ll ' 



where Cj is the coefficient defined in (18 1. We see that the direct (naive) evaluation of this 
expression at g = expiTr/ (£+ 2) is problematic as the some of the g-integers in the denominator 
become zero. Yet, writing out the products it is easy to show that for £ — 1 and £ = 2 this 
expression coincides with the modified local Hamiltonian h' defined in the last section provided 
that we set q — e"^/(^+^). For t — i, one checks that both Hamiltonians coincide at g = expi7r/5 
by using the explicit expressions given in |65j . In fact, we verified through explicit numerical 
comparison up to ^ = 6 the following statement: 

Conjecture 3.1 We have 

h' = hxxZe/2 g = exp {in/ (^ -f 2)) . 
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Wc have little doubt that this relation holds for general i. Hence, we propose the £-tli 
trigonometric model as lattice equivalent of the £-th model in the N — (2, 2) superconformal 
series with has central charge ci = Zljii + 2). This is compatible with the central charges for 
spin— £/2 XXZ chain [55J[B7]- The genuinely new feature of the present investigation is the 
construction of the two copies of the M = 2 supersymmetry algebra on the lattice. Moreover, 
the supersymmetry is not only a feature of the Hamiltonian: the supercharges have definite 
commutation relations with the transfer matrix of the fused vertex models as can quite directly 
be seen from the fusion equations and the algebraic Bethe ansatz |SH]. Eventually, notice that 
all the preceding expressions make formally sense in the limit £ — > oo. This is equivalent to the 
limit g — >■ 1 in the above expressions where [n] — >■ n. On the field-theory side the limit £ — >■ oo 
gives raise to a non-rational superconformal field theory studied in [SS] . It was shown in [70] 
that this type of non-rational limit can be related to free theories with continuous orbifolds. 
The fact that in the limit ^ — > oo we have a continuum of admissible twist angles, as can be 
seen from (14 1, supports our conjecture. 



3.5 Witten index 

We would like to understand if the models defined here above have zero-energy states. In general, 
if the Hilbert space % is finite, then a sufficient criterion for the existence of supersymmetry 
singlets is a non- vanishing Witten index |56j 

Indeed, as all excited states come in doublets with fermion number differing by one they cancel 
out of the trace. Only zero-energy states contribute. Therefore we find W — tr^(— 1)^ and thus 
\W\ provides a lower bound to their number. In the present case however, the dynamic nature 
of the supercharges makes a direct evaluation of the Witten index as defined above difficult 
because T-L in infinite-dimensional. Yet the trigonometric models allow to take advantage of the 
conservation of S\, and introduce and compute a Witten index- type object [I]. 

Recall that the trigonometric supercharges Qn at spin— £/2 change the total magnetisation 
by + 2)/2. Hence we conclude that the value of J = -I- 2)N/2 + Sf^ is conserved by 
their action. Moreover, we introduce the total particle number M = J — N (this is the actual 
particle number in the fermion model). Then for a doublet Qa? IV')) we find the values 
(M, M — 1). Hence we may introduce a version of the Witten index, given as the trace of 
(— 1)^^ over subspaces of H with constant J: 

oo oo 

Wj^Y. trw„(-l)*^<5M+iv,j = (-1)^ ^(-l)^tr„„,5M+Af,j. (25) 

N=l JV=1 

Thus, we have to evaluate the number vmn of states in Hn with total particle number M. 
We introduce its generating function fN{z) = X]m=o ^mnz'^^ with respect to M (for given TV 
the sum is actually finite). Assuming convergence, we can then write the generating function 
yV(z) = /Lo Wjz^ for the Witten index ( pS] ) in terms of the infinite series 

oo 

W{z) = J2 fA-z)z'. (26) 
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Open chains. The trigonometric models may be defined on an open interval because their 
supercharges are a solution of ^ with \x) = 0. For an open chain it is straightforward to 
evaluate the dimensions of the relevant subspaces as no translation invariance needs to be taken 
into account. The generating function is given by the simple expression 



fN{z) = 



1 - z 

Thus, we find the generating function for the Witten index to be 
wherefrom it follows that 

J = mod £+2, 
J = 1 mod £ + 2, 
otherwise. 




Closed chains: periodic boundary conditions. This case is more difficult as one has to 
take into account translation invariance. More precisely, we would like to count the number of 
states with momentum zero for N odd, and momentum tt for N even, with total particle number 
M fixed. Each such state can be generated from a representative = |mi, . . . , rriN) such that 
^f=i "^i ~ ^'-^^ ^ '-"^'^ j^^t have to count the number of inequivalent representatives 
where two configurations are considered to be equivalent when they are connected by a 
rotation. For iV even we have an additional constraint. In fact, each simple configuration 
l/Lt) = I mi, . . . ,TOAr) possesses a symmetry factor which is the smallest non-zero integer r such 
that T^\ijl) = l/x). The requirement to have states with momentum tt restricts r to be even. 
Therefore, we count the number of inequivalent representatives in this subset of configurations. 
This is a classical enumeration problem which can be solved by applying Burnside's lemma to 
the cyclic group. We defer the details to appendix [C] There, we derive the generating function 
for the dimensions vmn'- 

m=0 ^ 

Here gcd(a, b) denotes the greatest common divisor of a and b. For even N the factor 
= (— 1)™ subtracts the number of configurations with odd symmetry factors. 
Notice in particular that for z — 1 we obtain the total number of configurations without 
restrictions on the particle number: /Ar(l) = N^^ Em=o((-l)^^n^+ i))gcd(Ar,m) Qj^g^ j^^^^ 



we compute the generating function of the Witten index using (26 1, and obtain the surprisingly 
simple result 

Conjecture 3.2 The generating function for the Witten index Wj is given by 
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We verified the first few hundred terms in the Taylor expansion for ^ = 1, 2, 3, 4, but have not 
succeeded in finding a general proof. If this conjecture holds, we have 




0, J = mod £+2, 

(— 1)'^+-'^, otherwise. 



Both examples show that zero-energy ground states do exist in these models. In the case 
of periodic boundary conditions, it is easily verified that the result for the Witten index in 
compatible with the existence oi £ + 1 ground states at Sff — —£/2,—£/2+ 1, . . . ,£/2 — l,£/2 
if the number of sites is odd, and none if N is even. 



4 Spin— 1 models 

In this section we present different spin chains for spin one with lattice supersymmetry. These 
were obtained by solving (|6| with various additional requirements. The aim here is not to 
exhaust all possible solutions but rather to present some examples related to known spin chains. 
We start with the Fateev-Zamolodchikov spin chain and its elliptic generalisation in section [4?l] 
We point out that it possesses a supersymmetry any choice of its coupling constants. Moreover, 
we discuss some special points in the space of couplings were additional features emerge. In 
section |4.2| we proceed with the analysis of a supercharge which generates a spin chain related 
to the supersymmetric t—J model. In particular, we will present a supersymmetry-preserving 
elliptic extension of this model. In section [4. 3 1 we discuss a spin chain with particle number 
conservation mod 3 which is related to combinatorial problems, and interpolates between the 
trigonometric models for £ = 1 and £ — 2. 



4.1 The (elliptic) Fateev-Zamolodchikov chain 

In this section we investigate in detail the spin— 1 chain introduced by Fateev and Zamolodchikov 
[71] . and its elliptic extension [S3]. Heuristically this is motivated by the following observation. 
The trigonometric model at spin— 1 coincides with the twisted Fateev-Zamolodchikov chain at 
the special point 77 = 7r/4 as we saw in the last section. Its scaling limit is described by the 
second member of the JV ~ (2, 2) superconformal minimal model with central charge c = 3/2. 
This one is known to possess also an A/" = (1, 1) superconformal symmetry [75]. Indeed the 
relation between the Fateev-Zamolodchikov chain and supersymmetry in the scaling limit has 
been discussed in the literature [TT1[73]. We show that this symmetry exists on the lattice, too. 
It is genuinely different from the one analysed in the last section. 

We proceed as follows. First, we will prove the existence of the lattice supersymmetry in 
section [4.1.1[ and show that it is present even off the critical point. The sections |4.1.2| and 
|4.1.3| deal with some local and global symmetries of the Hamiltonian, and their interplay with 



the supercharges. The existence of supersymmetry singlets is discussed in 4.1.4 In section 
|4.1.5| we make a series of conjectures on the nature of these zero-energy states: in particular, 
we conjecture a relation to the weighted enumeration of alternating sign matrices. Eventually 



we come back to the point 77 = 7r/4 in section 4.1.6 where we show that the J\f = (2,2) 
supersymmetry is preserved off criticality. 
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4.1.1 Supercharges and Hamiltonian 



The type of Hamiltonian which we would like to study is locally of the form ( 19 1, and therefore 
HN^}2\j2ja + 2{s^r) - ^ I . (27) 

= 1 a,b=l 

Here Ja and Aat, a,h = 1,2,3 are constants to be fixed. The operators s", a = 1,2,3 are the 
spin-1 representation of su(2). In the usual basis where is diagonal they are given by 3 x 3 
matrices 





s 



= V2 




with = ± is^. Thus we have s^\m) = (to — 1)|to), and furthermore s^|0) = \/2|l) and 
= V2|2). 

We want to show that there is a local supercharge q which generates the local Hamiltonian 
density of (p7|. In order to construct it, we need the three states 



|Vo) = |10>-|01>, |V^i> = |20)-|02), |^2) = |21)-|12), 
which span a spin— 1 module of su(2) on two sites. The local operation q is defined through 
q|0) = xlV'o) - 2/|V'2), 

q|i> = IV'i), 

q|2) =a;|V^2)-2/|V'o). 

An explicit check shows that it leads to a nilpotent Qtv as ([6]) holds with vector |x) = 
(x^ - y2)|ll) - a;(|20) + |02)) - y{\22) + |00)). Writing out the Hamiltonian density h, we thus 
identify the couplings constants. We find that Aab — ^6a- The diagonal elements are given by 
A-aa = J a- The remaining couplings are quadratic polynomials in x^y: 

Ji = 1 - 4x?/, Jz = 1 + ^xy, J3 = 2{x^ + y^) - 1, 

Ai2 = l-^V\ ^i3 = (2x-l)(l-22/), A23-(2a;-l)(l + 2y). 

The Hamiltonian density h is invariant under spin reversal as it should be because (R® K)qR = 
— q. This implies that no second copy of the supersymmetry algebra can be constructed from q 
through spin reversal. From (|3| we conclude that the only possible twist angles are = 0, tt 
(the latter will play an important role in what follows). The case y = 0, which corresponds 
to the trigonometric limit, however admits arbitrary twist angles, and the supersymmetry 
is realised provided that ^ holds. In terms of the spin chain, twisted boundary conditions 
correspond to 

The fact that is arbitrary comes from the property that the local supercharge q is neutral for 
y = 0: it does not change the magnetisation unlike the trigonometric supercharge at spin one 
defined section [3l 
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Elliptic parametrisation, integrability. Fateev showed that the Hamiltonian (27) is 
integrable on along a two-parameter subnianifold in the space of couphngs |63j . The coupling 
constants are parametrised by Jacobian elliptic functions. The Ja take the values 

2 2dn27 _ 2cn27 

J^ — ; — Ji — ; — J?, — 



l + dn27' l+dn27' l+dn27 

Furthermore, the symmetric 3x3 matrix A has diagonal elements Aaa = J a, and off-diagonal 
elements 

2sn27 /cn7 cn27\ 2sn2^ f \ 1 

^12 = — , — ^ 1 : — ^ I , ^13 



A 
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1 -f dn 27 \ sn 7 sn 27 / 1 dn 27 \ sn 7 sn 27 

2 sn 27 / dn 7 dn 27 



1 -|- dn 27 \ sn 7 sn 27 



Here the two relevant parameters are thus 7 = 2_ftr?7/7r, where < 77 < tt is the so-called crossing 
parameter, and the elliptic name p = exp^—nK'/K), where K,\K' are the quarter-periods of 
the elliptic functions [74 . We call the generic case with non- vanishing nome q elliptic. The 
trigonometric limit corresponds to p — > (where K — > 7r/2): the coupling constants reduce to 



the values ( 20 ) 



In fact, it is just a matter of comparison of coefficients that the choice 
(1 -t- dn7) sn27 (1— dn7)sn27 



X 



2(1 -f dn27)sn7' ^ 2(l+dn27)sn7 



leads to supercharges that generate the Hamiltonian density of the spin chain ( 27 1 along the 
integrable line. Hence we conclude that the Hamiltonian of the (elliptic) Fateev-Zamolodchikov 
spin chain is part of an A/" = 2 lattice supersymmetry algebra for any choice of the parameters x 
and y, in sharp contradistinction to the XYZ chain at spin— 1/2 where the lattice supersymmetry 
is only present at a the special point r] — tt/3. 

It is interesting to compare this observation to the field-theory limit of the chains. In the 
trigonometric case, it is known that the low-energy states in the antiferromagnetic regime (the 
one studied here) are related to the spectrum and field content of an A/" = (1,1) superconformal 
field theory with c = 3/2 In the scaling window around the trigonometric point the 

appropriate off-critical extension thereof is given by the so-called super sine-Gordon field theory. 
The field theories possess an A/" = 1 supersymmetry for both left-movers and right-movers - a 
concept which does not exist on the lattice. We have therefore a lattice Af = 2 symmetry as 
predecessor of the field theory symmetries. 



4.1.2 Rotational invariance and duality transformations 

Motivated by the discussion on duality transformations for the spin- 1/2 XYZ chain at the 
combinatorial point given in section |2.4[ we discuss the symmetries of the present supercharge 
with respect to rotations. As before, we denote by p'^{9) = exp (i0s°) a local rotation about 
the a— axis by an angle 9, where s° is one of the su(2)-generators in the spin— 1 representation, 
and by ^%[{0) — IljLi Pj{9) a global rotation. We start with generic values for x,y for which 
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no simple rotational invariance is expected. Writing q = q{x,y) wc find that 

(pi(7r/2) ® pi(^/2)) q{x, y)p'{-n/2) = {x - y)q |(|^^) > (28a) 

(p2(vr/2) ® p2(^/2)) q(x, y)p\-7:/2) = [x + y)q lo^T^) ' ^^^^^ 

(p3(7r/2) ® p3(^/2)) q(a;, y)p3(-7r/2) - q (x, -y) . (28c) 

In the case of untwisted boundary conditions = we conclude thus that the Hamiltonian 
transforms according to 



n],{^/2)HN{x, y)n],{~^/2) ^{x- yfH^ 

nl{^/2)HN{x,y)nl{-n/2) = {x + yfHN 
n%{7r/2)HN{x, y)n%{-7r/2) = Hn {x, -y) . 



y + I X + y — I 



2{x^y) 2{x-y) 
X — y + \ X — y — 1 
2{x + y) ' 2{x + y) 



We see that up to rescaling the spectrum is invariant under the transformations in parameter 
space induced by these rotations. Unfortunately, if we choose the twist angle 4> = n, which will 
prove to be the most interesting choice, then only the third transformation survives. We discuss 
briefly some special points: let us suppose that x is left invariant by (28a I. Then we have either 
(i) X = y+1, which can be transformed by (28b I according to 1, y) (l/(2j/+ 1), 0) or (ii) 
X = —1/2. If we set 2/ = C/2 in the second case then we see that (28a I leads to ( (3— C)/(1 + C)i 
and (28b) to C — > (C + 3)/(C — !)• We recognise here the homographic transformation which 
appears in the spin— 1/2 XYZ chain at the combinatorial point as discussed in section [2!4| 

Eventually, there is a symmetry relating x and ~x. Consider the projector = s^(l =ps'^)/2 
on the spin states 2 and respectively. It follows that 

(e-P* ® e-P*) q(x, y)e"'^P^ = q(-x, y), 

what implies that the Hamiltonian enjoys the transformation property 

N 

n%HN{x,y){n%)-' = HN{~x,y), = l[e'^pf. 

All these symmetries are simple in the sense that they can be built as tensor products of 
local transformations. Less obvious seems the following statement: 

Conjecture 4.1 For the twist angle (/) = tt the spectrum is invariant under the exchange of x 
and y. 

This statement implies the existence of a unitary transformation V(x,y) with the intertwining 
property V{x, y)Hi^{x, y) = His[{y, x)V{x, y) with V{x, y)~^ = V{x, j/)* = V{y, x). The explicit 
construction of the simplest examples with N — 2 and = 3 sites shows that V{x, y) is not a 
tensor product of local transformations, even though the eigenvalues of the Hamiltonian density 
are invariant under the exchange of x and y. 
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Finally, we turn to the trigonometric limit y — 0. We have a full invariance with respect 
to rotations about the z— axis: {p^{9) (8) P^{0)) q P^{—(^) = q- This is nothing but the U(l)- 
symmetry which leads to the conservation law for the total magnetisation in the Fateev- 
Zamolodchikov chain. It becomes a full SU(2) symmetry at the point x — l,y — which 
corresponds to the XXX chain at spin— 1 [751 [7B]. Indeed, the supercharges itself are invariant 
because (g) p''(6')) q(l, 0) p"-{-0) = q(l, 0) for all a = 1, 2, 3 at this point. The invariance 

is due to the fact that the supercharge maps a spin— 1 module onto a spin— 1 module, the latter 
being realised on two sites. 



4.1.3 Local quantum group symmetry 

The fact that the local supercharge q makes explicit the SU(2) symmetry at a; = 1, y = leads 
quite naturally to the question if for x — {q + q^^)/2, y — the local quantum group invariance 
of the Fateev-Zamolodchikov chain may be seen from the supersymmetry. This is indeed the 
case as we show now. To this end, we consider the gauge-equivalent local supercharge 

q' = q+q0, q^|m) = |/3) ® |m) + |m) ® |/3) with \13) ^ ^{q-^ ~ q)\l) . 

As noted before, we add only a redundant part which has no global effect. The new local 
supercharge acts according to 

q'|0>=<7-i|10)-g|01), 

q'|l) = |20)-|02) + (g-i-g)|ll), 

q'|2>=g-i|21)-g|12). 

The new supercharge q' solves ^ with \x) — ~{q\02) + q~^\20) — |11)) which is the quantum- 
group singlet. Moreover, on the right-hand side of the preceding equations, we recognise 
basis states for the spin— 1 module of Uq(sl2) realised on two sites. We thus see that here 
the supersymmetry is a consequence of the structure of the g-Clebsch-Gordon coefficients. 
Moreover, the given structure implies immediately that q' commutes with its generators in the 
sense that 

q's± = A(s±)q', qV = A(s3)q', 



where A is the comultiplication defined in ( 22 1 . According to ( 10 1 the corresponding Hamiltonian 



densities h and h' differ by a local boundary term. Indeed, we find 

— 2 2 

h' ^ (s^ ® 1 - 1 (g) S^) . 

This is precisely the Hamiltonian density of the quantum group invariant Fateev-Zamolodchikov 
spin chain studied in |77j . 



4.1.4 Supersymmetry singlets 

Given that the elliptic Fateev-Zamolodchikov spin chain has a supersymmetry for any choice 
of the parameters x,y we ask if any solutions to the equations Qn\^) = and (3]y_-^|<E>) = 
exist, and how the space of solutions depends on the parameters. We use the conjugation 
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technique described in section |2.4| in order to relate different values of cc, y. Consider the 
operator m — m(A, 9) defined through 

m(A,^?)|0) = (cosh (^0 |0) +sinh (^0 |2) 
m(A,0)|l) = |l), 

m(A,^?)|2) = (^sinh Q |0) +cosh (^0 |2)) . 

It follows that 

(m(A, 6) ® to(A, 61)) q(x, y)m{X, Oy^ = A~iq(a;', y'): 

where 

We thus see that A plays the role of a dilatation parameter, whereas 9 is the angle of a "Lorentz 
transformation" on the parameter space. We allow negative values for A, and thus include 
reflections with respect to the origin. As in the case of the XYZ chain studied in section |2.4[ 
the transformation can be promoted to the full chain: we have Mjv(A,6') = m{X,9)^^ . The 
conjugation of the supercharges gives 

MN+iiX,9)QN{x,y)MN{X,9)-' = X-^Qn{x' ,v'). 

Notice that Mn{X,9) preserves the momentum spaces where TM{(t>) = (—1)^+^ for = 0,7r. 



Hence, from the conjugation argument given in section 2.4 we conclude that Hj\;(x,y) and 
i?Ar(x', y') have the same number of zero-energy states. Thus, we need to understand which 
points can be mapped onto each other. It is easy to see that (a;')^ — (y')^ — X^{x^ — y^). This 
divides the (x,y)-plane in (i) time- like points x'^ — y^ < 0, (ii) light-like points x^ — y^ and 
(Hi) space-like points x^ — y^ > with respect to the origin x — y ^ 0. The latter plays a 



distinguished rule, as it is the only fixed point. The transformation (29 1 preserves the nature 
of the points, i.e. if {x,y) is time-like then {x',y') is, and so on. It follows that for all time-like 
choices of the coupling constants x, y, the Hamiltonian has the same number of zero-energy 
states. The same applies to all space-like choices for x,y. The conjecture |4. 1| implies then that 
for all choices with |a;| ^ \y\ the Hamiltonian has the same number of supersymmetry singlets. 
Therefore, it is sufficient to analyse a single, not light-like point in the (x, y)-plane. As opposed 
to the XYZ chain studied in section [2^ there does not seem to be a simple choice however. We 
content ourselves with the following 

Conjecture 4.2 For any number of sites N and twist angle cf) — n there exists a choice of 
parameters x,y with x ^ y so that the Hamiltonian has a single zero-energy state in the sector 
Sfj = mod 2. // y = the single zero-energy state has Sf^ — 0. Furthermore, for the twist 
angle = there is a choice of parameters x, y with x ^ y so that the Hamiltonian has no 
zero-energy states. 

If this conjecture holds then it implies the existence of exactly one zero-energy state for all 
choices with |a;| ^ \y\ if the twist angle is set to = tt, and none if </> = 0. The existence of a 
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singlet for x 0,y — for the twisted chain was already observed in [78, where a staggered 
spin chain based on the quantum group Uq(s[(2|l)) was studied. Its spectrum contains the 
spectrum of the twisted Fateev-Zamolodchikov spin chain. 

Next, let us consider the light-like points with |x| = |?;|, x ^ 0. We show in this case 
that there is at least one ground state. We start with even N. Invoking the conjecture |4.1| 
we see that it is sufficient to consider the half-line x — y, x > 0. All choices for a; > lead 
to the same number of singlets. The point x — y — 1/2 is particularly convenient because 
most of the terms in the Hamiltonian cancel as can be see from section |4.1.1[ From any spin 
configuration n which contains only O's and 2's we can build a = 12j=o^ ^atI/^)' which 
is an eigenstate T;v((/) = tt) with eigenvalue (—1)^+^ = 1. If restricted to this subspace, the 
Hamiltonian becomes very simple. Along each row and column, the sum of its entries are zero. 
Hence, up to a multiplicative factor it has the zero-energy state |^) = \fi) where fi runs 
over all spin configurations which do not contain any I's. Next, let us turn to odd iV = 2n -|- 1. 
In this case it turns out to be convenient to choose the half-line x = —y, x > 0, and pick 
the point x — —y — 1/4. Instead of we consider the unitarily equivalent Hamiltonian 

H'jq = UHnU^^ with U = VW where W = exp (^^t^ J2f=i J ■^'j^ ■ ^ is a diagonal operator 
defined as follows: denote by the number of of I's in the spin configuration /i, for 

example i/i(0100221) = 2, then V\^) = i'^^'^^'^l/x). is invariant under translation, and we are 
interested in its zero-energy states in the translationally invariant sector. It turns out that if 
we set X = 1/4, y = — 1/4 then the corresponding Hamiltonian density has again the property 
that the sum of its components along each row or column is zero. This time, no restriction to 
any other subspace than the translation sector is necessary. Hence, a zero-energy state of H'j^ 
is given by |^') = what implies that Hn has at least one zero-energy state. Using the 

conjugation argument we conclude that for all choices of x,y with \x\ — \y\ and x,y ^ there 
is at least one zero-energy states in the momentum sectors of interest. 

We are now left with the point x = y = 0. It is known that in this case the Fateev- 
Zamolodchikov Hamiltonian can be mapped onto the supersymmetric t — J model. This case 
will be analysed in section (4.2 1, where we will show that Hm possesses zero-energy states 
even for = 0, and thus demonstrate that their number may change when parameters of the 
Hamiltonian are varied. 



4.1.5 Combinatorial nature of the singlet states 

It is natural to ask if for the twist angle (j) — tt the singlet state and its components, defined 
through 1^*) — where /i runs over all spin configurations, possess interesting properties. 

Choosing its components to be coprime polynomials in x, y, an example of a remarkable sum 
rule is the following observation: 

Conjecture 4.3 The square norm of the supersymmetry singlet eigenvector factorises into 
two polynomials in x,y. At at odd N, one of these polynomials is the component ipii...i. 



Trigonometric case. We observed a refinement of this conjecture in the trigonometric case 
when y ~ from exact diagonalisation. To state it, we need to fix the normalisation of the 
ground state vector j^"). As the Hamiltonian is a quadratic polynomial in x we may choose its 
components to be coprime polynomials in x. This fixes the state up to an overall multiplicative 
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constant. We fix it through the asymptotic the hmit x oo. The Hamiltonian reduces to 



N 

lim x--'H^ = 2 ^(44+1 + {s'^r + (sUi)") (30) 



.7 = 1 



and therefore does not depend on the twist. It is easy to see that the Hamiltonian is positive 
and annihilates the states |11 • • • 1) for odd N, and |0202 • • • 02) + |2020 • • • 20) for even N. These 
are the only zero-energy states on the sectors where T^icf) = vr) = (—1)^+^. It follows that for 
N odd the component V'li - i must have the highest degree in x amongst all the components of 
the singlet vector if they are chosen to be polynomials. The same applies to the components 
V'0202 - 02 and '02O2O- -2O for even N. We fix the normalisation such that 



V-ii.-i (2x)''"(l + o(l)) for TV odd, (31a) 

a;— >oo 

V'0202...02 - (2a;)'*"(l + o(l)) for iV even. (31b) 



Conjecture 4.4 (Degree) With this choice of normalisation all components are polynomials 
in X with integer coefficients if we adjust the degree to 




71^, 7V = 2n+l, 

71(71-1), N^2n. 



We are ready to analyse the states in a systematic way. It turns out that both their 
components, and square norms are related to problems of enumerative combinatorics. This is 
rather similar to the case of the XXZ chain with anisotropy A = —1/2 whose combinatorial 
features were first observed in [301 ST] . In order to understand the combinatorial nature of 
the supersymmetry singlet in the present context, we need some basic notions on the weighted 
enumeration of alternating sign matrices. Consider the set of 71 x n alternating sign matrices. 
Their entries take the values or ±1. The first non-zero element in each row and column is 1, 
and the non-zero elements alternate in sign. We assign the weight t'' to a matrix if it contains k 
—Vs. The generating function for given size is A„{t). The example of 3 x 3 matrices is shown in 
figure [2] and gives A^lt) ~ 6 + t. Next, restrict to (2n 4- 1) x (2n + 1) alternating sign matrices 

/ioo\ /ioo\ /oio\ /oio\ /ooi\ /ooi\ /o 1 o\ 

010 001 100 001 100 010 1-11 
Vooi/ Voio/ Vooi/ Vioo/ Voio/ Vioo/ Vo 1 0/ 

1-fl + l + l + l + l+t = 6 + t 



Figure 2: All 3x3 alternating sign matrices: six permutation matrices with weight 1 and a 
single matrix with weight t. The generating function for these is ^3(t) = 6 + t. The last matrix 
is the only 3x3 vertically symmetric alternating sign matrix. Therefore we have Ti(t) ~ 1. 

which are symmetric about the vertical axis. Assign weight t'^ to any such matrix if it contains 
k — I's within the first n columns. Then, the generating function is given by Tn{t). From figure 
[2] we infer that for n = 1 there is only one vertically-symmetric alternating sign matrix which 
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gives Ti(t) = 1. In general, both T„(i) and A„(t) are given in terms of determinantal formulae. 
The relevant polynomials can be found in Robbins' [7S] and Kuperberg's work [55]. We need 
the following expressions: 



2n-2 



o<i.j<n ^—^ \k — il \2j — k 

2n~l 

Rnit)^ det y r(i,fc,i)y(j,fc,o)<2j+i-'=. 

0<i,?<n ^ — ^ 
fc=0 

The coefficients Y{i,j,k) are 

i + k \ / i + k+1 
+ 1 + fc - j7 ^ V2i + 1 + fc - J 

The polynomial An{t) is then given by 



Yit,j,k) 



A„(t) 



T„{t)R„{t), N^2n+1, 



These polynomials appear in the supersymmetry singlet provided that we identify t = Ax^. 
Indeed, up to iV = 11 sites we checked the following statement: 



Conjecture 4.5 Choose the normalisation as in (311. For odd N — 2n+l we have 

whereas for even N — 2n 

V'0202...02 = r„(4a;2), \\^\\^ = A2n{ix^). 

We see that in both cases there are distinguished components related to vertically-symmetric 
alternating sign matrices whereas the square norm gives the polynomials for unrestricted 
alternating sign matrix enumeration. For odd 2n + 1 we found another sum rule related 
which is related to the polynomials To state it, we transform the state according to 

In fact, U is the unitary transformation which spreads out the twist over the whole system so 
that for N odd UHnU^^ is invariant under translations. Next, we denote by i'o(m) the number 
of spins with value in the configuration /i. For example fo(0100221) = 3. If the number of 
sites is odd = 2n + 1 and the states under consideration have 5*^ — the 1 < foifi) < n. 
With these preliminary definitions we claim the following: 

Conjecture 4.6 For the components of the transformed states defined through |$) = 
we have the linear sum rule 
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Elliptic case, odd length. It is quite natural to ask how the combinatorial content gets 
deformed as we go off the critical line y — 0- We have not been able to formulate a conjecture 
on this deformation except for the point r] — 7r/3, corresponding to x — —1/2 and y = (/2 
which we encountered already when studying the duality transformations. At this point we 
find the coupling constants: 

Jl = 1 + C, ^2 = 1 - C, ^3 = Ai2 = -2J3, Ai3 = -2J2, A23 = -2Ji. 



As in the spin-1/2 case (12 1 the Hamiltonian is a quadratic polynomial in (. Thus the 
components of the singlet for (j) — n can be chosen as mutually coprime polynomials in C. This 
fixes the normalisation up to an overall factor. We investigated this state up to = 7 sites, 
and put forward the following 

Conjecture 4.7 For rj = 7r/3 and N — 2n + 1 sites we find in suitable normalisation the 
component: 

^ii...i(C)=(^— j P—^[sTc 

where the Pniz) are related to the polynomials found by Bazhanov and Mangazeev \5(A \51'l for 
the spin— 1/2 XYZ chain (we recall them in appendix\B\ see also f5S^ ). The second factor in 
the square norm can also be written in terms of the polynomials PniC,) '■ 



3 + C 

For N — 2n sites, the normalisation can be adjusted in such a way that the square norm is 
given by 



4.1.6 Extended supersymmetry 

One of the main motivation to study in detail the Fateev-Zamolodchikov spin chain was the 
coincidence of the second member oi Af — (2,2) superconformal unitary series with a,J\f= (1,1) 
superconformal field theory. In terms of the lattice models, we understand now that in the 
trigonometric case the supersymmetry is enlarged if we set 77 = 7r/4. In this section we show 
that this feature survives off the critical point. As a heuristic motivation, let us recall some 
results on the super-sine-Gordon theory which describes the scaling limit of the spin chain. It 
was found in [501 EI] that the soliton scattering matrix of this field theory factorises according 
to S{9) = iS^q'q (6*) (8)5sg(^: /3sg) where 9 is the standard rapidity variable. The factor S^gQ{9) 
is the iS-matrix of the so-called K-th restricted sine-Gordon field theory, and here only relevant 
insofar as for K = 2 it corresponds to the Af = (1,1) supersymmetry. The second factor 
Ssg{&, Psg) is the scattering matrix of the usual sine-Gordon field theory. The latter is known 
to admit anAf= (2, 2) supersymmetry if the sine-Gordon coupling takes the value /3gQ = 167r/3, 
as was pointed out in |39j . This point is equivalent to 77 = 7r/4 in the lattice model. 
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Guided by this heuristic consideration we seek for an elHptic extension of the spin— 1 
lattice Af — (2, 2) supersymmetry algebra presented in section [s] There, we found two 
local supercharges q and q which changed the magnetisation by AS^ = — 1 and AS^ = +1 
respectively. Changing this to AS^ = ±1 mod 2 = 1 mod 2 destroys unfortunately the natural 
distinction of q and q (this is a general phenomenon for even £). We circumvent this problem 
by constructing local supercharges with built-in spin-reversal symmetry {R ^ R)q^ R — ±q^ 
with the requirement that they both generate the same local Hamiltonian — . 

We start with a supercharge with the property (i?® R)q^R — — q^: we determine constants 
ai , a2 , 03 , 04 such that the action 

q-|0) =a3(|21) + |12)) -02(101) + |10)), 
q-|l>=ai(|22>-|00» + a4(|20>-|02)), 
q-|2) =02(121) + |12))-a3(|01) + |10)) 

leads to a nilpotent Qn- We prefer to keep 02, 03 7^ in order to assure that in the trigonometric 
limit the solution reduces to a suitable odd combination of the trigonometric supercharges. 
Now, one has to construct a vector \x) such that (|6| is solved. This is indeed possible provided 
that 04 = 0. Thus, we have three free parameters oi, 02, 03. 

Next, let us construct a supercharge such that (i? (8) i?)q+i? = q^. In components, we write 

q+|0> =63(121) + |12)) + 62(101) + |10)), 

q+|l) = 6i(|22) + |00)) + 64(120) + |02)) + 65III), 

q+|2) =62(121) + |12))+ 63(101) + |10)). 



In fact, it contains a redundant part in the sense of proposition 2.1 which can be eliminated 
through a redefinition 62 — >■ 62 -I- 65. Hence we may safely set 65 = 0. The requirement that this 
local operation leads again to a nilpotent supercharge leaves us with the additional conditions 
6162 = 6364. Hence, we are left with three free parameters 61, 62, 63. 

Eventually, we demand the Hamiltonian densities associated to the two supercharges coincide: 
/i+ = . The system of equations for the remaining free coefficients Oj , 6j, j = 1, 2, 3, can be 
solved by 

/— Q'lQ'2 ^1 2 r»/ 2 2\ 

61 = V2a3, 62=-^, 63 = -^, 01 = 2(02-1-03). 

As we are free to rescale the parameters we set oi = 1. We parametrise the remaining constants 
by a parameter — 1 < C ^ 1 such that 



= — ^= , 03 = 



2^/2 ' ' 2V2 
This fixes all coefficients as functions of Moreover, the q^ lead to supercharges with the 



mutual anticommutation relations of the J\f = (2,2) supersymmetry algebra (111. As the 
particle number conservation holds for generic values of C only mod 2, the only possible twist 
angles are cj)^ = 0,7r. The local Hamiltonian /i+ = generated from these two supercharges 



corresponds to a particular version of the local Hamiltonian h oi ( 27 1 up to a similarity 
transformation 



/i±= 1 
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provided that we set 



1 - dn 27 _ 2Kri _ K 
^ ^ 1 + dn 27 ' 

This result is the elHptic version of (21 ): the fact that the both spaces acts on are transformed 
differently implies that is a twisted version of Hj^i. With the unitary transformation 

U — exp (^iTT js^) we find = UHnU~^. We denote by (j)± the twist angle for the 
present supercharge, and by </> the twist angle for the supercharge of (27 1. Then we obtain that 
0± = + Nt: mod 27r. 

We have thus established that at 77 = 7r/4 the twisted Fateev-Zamolodchikov spin chain 
possesses an A/" = (2, 2) lattice supersymmetry even away from the critical point. Given this 
fact, we may ask if it leads additional zero-energy states present at this point. The exact 
diagonalisation for = 3, 5, 7 sites suggests that the answer is affirmative: 

Conjecture 4.8 For an odd number of sites N = 2ri + 1, twist angle (j) = n and crossing 
parameter 77 = 7r/4 the Hamiltonian possesses two additional zero- energy states in the sector 
with Sf^ — 1 mod 2. In the trigonometric limit these states can be found in the sector where 

s% = ±i. 

Notice that this is not in contradiction with the general discussion in section |4.1.4| whose 
considerations applied only to the sector with Sf^ = mod 2. 



4.2 Spin chains and the supersymmetric t — J model 

For general spin— £/2 there is a simple case where ([6| is satisfied with \x) =0: suppose that q 
is only non-zero when acting on a fixed component |m). Moreover assume that the resulting 
states \ jk) are all such that j, k ^ m. One encounters such a case for the Fateev-Zamolodchikov 
chain at the point 77 = 7r/2, i.e. x = y ^ Q, were the supercharge acts like 

q|0)=0, q|l) = |20)-|02), q|2) = 0. 

The local Hamiltonian is then given by /i = 1 + 11 where 11 is the graded permutation operator 
defined through 

n|l, 1) = |1, 1), n|m, 1) |l,m), n|m,m') = -|rn',m), to,to' = 0,2. 

If we consider the 0, 2 degrees as particles, then their individual numbers and n2 is are 
conserved. The local Hamiltonian permutes them like fcrmions with strong repulsion (i.e. a 
site can either be empty or occupied by a single particle) on a neutral (bosonic) background 
of I's. This is precisely the local operation of the supersymmetric t — J model [52] [S3]. The 
only difference here is that, as opposed to fermionic variables, the spin variables commute for 
different sites. In fact, the spin chain presented here differs from the supersymmetric t — J model 
by a generalised Jordan- Wigner transformation [84]. The dynamical lattice supersymmetry 
was found in |lj in the fermion language. Here we prefer to work in the spin language. 

We would like to study a simple deformation of this supercharge. Hence we choose 
q|0) = 0, q|2) = and q|l) = aoo|00) 4- a2o|20) -I- ao2|02) + a22|22). It is then trivial to see that 
^ holds. We restrict our considerations to a special case 

q|l) = |02)-|20) + C(|00)-|22)), ( G 
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We see that (i? (g) R)qR = — q and hence we have spin-reversal symmetry. We see that the 
Hamihonian generated from this supercharge is a quadratic polynomial in ^. Moreover, the 
equation ([3| shows that the only admissible twist angles for non-zero ^ are given by ^ = 0,7r. 
We focus on periodic boundary conditions and thus = 0. Switching on the parameter ^ 
breaks particle number conservation for the Os and 2s but not the Is. In general, we expect the 
singlet states to be rather complicated superpositions of polynomials in C (as is the case for 
the XYZ chain along its combinatorial line). Yet, we show here below that some of them are 
robust against this elliptic deformation. 



Witten index. We start the study of the ground state by evaluating a Witten index which 
is related to a conservation law like in the trigonometric models of section |3] Let us denote by 
Uj the total number of particles of type j = 0, 1, 2. The supercharge decreases ni — ni — 1, 
and increases the number of sites N ^ N + 1. Thus, J = N + ni is conserved, and we may 
introduce a Witten index Wj as in section |3.5| 



Wj^Y. ti-nA'^r'Sj^N+n, = (-1)^ ^(-l)^tr«„<5xw+„, 



N=0 N=0 



We denote again by h'ni,N the number of distinct momentum states with zero momentum as 
N is odd, and momentum tt for N even, containing ni particles of type 1. A slight variation 
of the argument given in section |3.5| and appendix [C] shows that the generating function 
fN{z) = J2'^,=o '^nuNz"^ is given by 



N-l 



ni—0 

Computing the generating function generating function >V(z) of the Witten index, we find 



again a surprisingly simple and regular series expansion as in section 3.5 Yet, as before we do 



not have a proof, but tested the first few hundred terms of the series expansion. We obtain 

oo oo 

W(z) = Wjz' = Y fj(-')^' = 

J=0 J=0 

and thus conclude that 

(32) 

There are thus at least two zero-energy states for given odd J = N + ni for periodic boundary 
conditions. 




Zero-energy states. As opposed to generic spin— 1 chains, the present model possesses zero- 
energy states which can easily be determined. We consider only periodic boundary conditions. 
We start with zero-momentum states at odd A''. Observe that any state which is only build 
from particles and 2 is annihilated by Qn- Therefore, we restrict our considerations to the 
subspace in which no local spin is in the state On this subspace the Hamiltonian density h 
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has two independent null vectors which are simple tensor- product states |x±) — \v±) l^i) 
with \v±) — (|0) ± \2))/y/2. There are thus two zero-energy states of Hn 



\'i'N) = \v±)^---\v±)- (33) 

These states are eigenstates of the spin-reversal operator: i?7v|^^) = and belong to 

the sector with 5^ = 1 mod 2. If C = then the Hamiltonian conserves Sf^. On the subspace 
considered here, we have Sff — —N, —N + 2, . . . ,N — 2,N . By projection of |^^) on these 
subsectors we obtain therefore iV -I- 1 zero-energy states for C = (not 2(A^ + 1) because of the 
spin-reversal symmetry) . 

Next, let us consider states with momentum tt at even N . We develop a suitable modification 
of the states for zero momentum by allowing a single site with state |1). Hence, no pairs 
|1) ® |1) are present. The Hamiltonian density at C = has the zero-energy states \v±), and 
|wo) = |01) — |10), \v2) = |21) — |12). Using this fact, it is easy to see that the states 

N 

\'^%) = Y^{-^y\v±)®---\v±)® ii) ®\v±)®---®\v±) 



are annihilated by the Hamiltonian with C = 0. Upon projection we conclude that there is a 
single momentum state in each subsector with = — A^-f-l,— Af + 3,...,A^ — 3,A^— 1, thus in 
total N zero-energy states for C = 0. If however, C 7^ this construction is not possible as the 
spectrum of the Hamiltonian density changes. Nonetheless, both cases C = and C 7^ lead to 
the same Witten index. To see this consider first C = and fix J = -f ni to be odd. We thus 
have J +\ zero-energy states for ni = 0, and J — 1 zero-energy states for ni — 1. This gives 
Wj = {J + 1) — {J — \) = 2. Next, fix J even, then Wj = trivially. For C 7^ there are only 
two zero-energy states at N odd with ni = 0. Thus Wj is the same as in the case of C = 0. 
Despite the simplicity of this example, it has an interesting feature: the zero-energy 



states (33) are present for any C, and thus stable against a certain supersymmetry-preserving 



perturbation away from the critical point. 



Conserved quantities. The supersymmetric t — J model is known to be integrable. In 
particular, its conserved quantities were systematically determined by Essler and Korepin [53] . 
It is thus natural to investigate the question if the integrability of the present model survives if 
the parameter ( is switched on. 

We give an argument in favour of this conjecture using Reshetikhin's criterion [5S] for the 
existence of a higher conserved quantity H^"^ — I^jj+i' ^j+ij+z]- Obviously, it commutes 

with the translation operator i?^'* — Tjq. The criterion states that commutation with the 
Hamiltonian Hlj' — is assured if there is an operator X acting on V (^V such that 

The right-hand side is a local boundary term in the sense defined above, and thus cancels out 
in the summation over all sites. In the present case, a non-trivial X can indeed be found: 

X = -2{1-C^)^{h + X'), 
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where X' itself commutes with h. In terms of the su(2) generators at spin-1 we find 



X' = (s^ «1 s" + S" (g) s+ + {s^ (g) S^, S+ (g) S" + S" «) s+}) . 

(1 - C^) ' 

It is known that the existence of if^"* imphes that there is another operator — [_B,i?^''], 

where B — X^jLiJ^jj+i is the boost operator. It commutes with 7?^' |HB] what follows 
from the Jacobi identity. The explicit construction of higher conserved charges is a non- 
trivial problem. Yet the simple result presented here makes plausible their existence. This 
problem might find an elegant solution by relating the present Hamiltonian to a transfer matrix 
constructed from an i?-matrix which will most likely be an elliptic extension of the solution to 
the Yang-Baxter equation presented in [85] . 

4.3 The mod— 3 chain 

In this section, we introduce another set of supercharges for spin one. So far, our discussion 
focused on trigonometric and elliptic supercharges. The latter conserve the number of particles 
mod 2. If ^ -|- 1 = mod p then it is possible to modify this to a mod p-conservation. Here we 
concentrate on £ = 2, p — 3. 

We present a solution of (|6]) with the requirement m -f 1 = j + fc mod 3 (and the condition 
that the components which violate exact equality are non-zero). It is parametrised by two 
coupling constants A, and its action is given by 

q|0) =C|22), 

q|l) = AC(|02) + |20))-|ll), (34) 
q|2) = A|00). 

Indeed, one may verify that the condition ([6]) holds if we set |x) — AC(|02) + |20)). Writing 
q — q(A,C) we have the property (i? g) i?)q(A, C)^ — q(Ci'^) under spin reversal. Hence for 
periodic boundary conditions the line A = defines a spin-reversal symmetric Hamiltonian. 
Moreover, we have q°P = q and therefore refiection symmetry of the Hamiltonian. The mod— 3 
conservation implies that the supersymmetry admits the twist angles — 0,±27r/3. Here we 
focus only on periodic boundary conditions. 

Asymptotic limits. Let us investigate some special limits in order to understand the 
structure of the ground states of the associated spin chain for an odd number of sites. As we 
may exchange A and C by spin-reversal, we focus mainly on C. We consider three cases: 

• At C = we find q|l) — — 111) and q|0) = 0,q|2) — A|00). Hence the supercharge acts 
trivially on and decouples it from |0), |2). Its action on the subspace spanned by the 
latter coincides with the one for the XXZ chain at A = —1/2. For N odd, there is a 
trivial ground state \^%s) = |11 • • • 1), and the two ground states of the XXZ chain, 1^*^), 
in the sectors with 5'^ — ±1, which can be mapped onto each other through spin reversal. 
These are the Razumov-Stroganov states [40j which can be normalised in such a way that 
all components are integers with the smallest one being 1, and the largest one An for 
A = 2rt + 1 where An is the number of n x n alternating sign matrices. 
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• If C — > oo wc find, upon appropriate rescaling, tfiat q|0) ~ |22), q|l) = A(|02) + |20)), and 
q|2) = 0. Up to a cyclic shift, given by the unitary transformation 




it coincides with the trigonometric spin— 1 supercharge of section [3] before imposing the 
spin-reversal symmetry. If A = l/v^ we recover the twisted Fateev-Zamolodchikov spin 
chain at 77 = 7r/4. In our previous discussion we argued that for TV odd there are three 
zero-energy states, one in each of the sectors S% = 0,±1. Hence, we see that in some 
sense the present model provides a supersymmetry-preserving interpolation between the 
first two members of the trigonometric models studied in section [3] 

If A = C = 1 then the model has a full cyclic symmetry: the Hamiltonian is invariant with 
respect to the cyclic shift (fi (g) f2)/i(f2~^ (g) Jl"^) = h. In fact, the Hamiltonian density 
becomes h = 2 — p with the operator p acting according to 

p\j,k)^\j + l,k-l) + \j^l,k + l), 

where we identify values of j, k mod 3. The Hamiltonian is trivial as it can be diagonalised 
by a simple Fourier transformation. We introduce the states 



2 iTr/6 



and build the Hilbert space for a chain with N sites from tensor products |j5i, . . . ,pAr). 
Then, the Hamiltonian = X]j=i ^j-j+i diagonal: 

The zero-energy states are thus given by the states with pj = p = 0, 1, 2 for all j = 1, . . . , A^. 
For iV odd, they are precisely the supersymmetry singlets (as the latter exists only in the 
momentum space with i^v — ("l)^"*"^)- Thus we find again three supersymmetry singlets 
for N odd, and none for N even. 



Zero-energy states. The last case in the preceding list is interesting because the general 



case can be related to it by conjugation, following the strategy outlined in section 2.4 We 
restrict our considerations to C, M > (other cases can be related to this one through simple 
unitary transformations). We introduce the diagonal matrix m = diag((A'^C)~^^^i li (AC^)~^^'^)i 
and find 

(m (g) m)q(A, C)m^^ = q(l, 1). 

Thus, using Mn = m®^ , one has Mn+iQn{\OM]^'^ = Qn{\ = 1,C = !)■ As we know the 
ground state structure in the case A = 1 , C = 1 in the supersymmetry-relevant momentum 
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sectors, we conclude that for odd N the Hamiltonian Hn generated from the supercharge 



(34) has three zero-energy states which are invariant under translation. They can be chosen 



according to the values of 5*^ = 0,1,2 mod 3. 

It would certainly be interesting to understand if the components of these vectors relate 
in any respect to enumeration problems. This is the case if C = (or A = 0) because the 
XXZ ground state at A = —1/2 relates to the alternating sign matrix numbers. For general 
values of A, C we normalise the singlets so that their components are coprime polynomials in the 
parameter. The inspection of small systems suggests that a combinatorial interpretation might 
not be straightforward, since not all the coefficients of these polynomials are positive (what 
is already seen for = 5 sites). Nevertheless, there seems to be an interesting underlying 
structure. The following conjecture, again based on exact diagonalisation up to = 7 sites, 
remind us of the sum rules for the Fateev-Zamolodchikov singlet: 

Conjecture 4.9 (Sum rules) With the given normalisation the square norm of the singlet 
state for odd N factorises into two polynomials in C, A with positive integer coefficients. Along 
the line X = we find more explicitly 



where vi{y) is the number of 1 's in the spin configuration n For the state we find along the 
line X = C the sum rule: 

^r(^)<^3^"c'"+v?i...i. 



5 Conclusion 

In this article we presented a general framework for dynamical supersymmetry of spin chains, 
and worked out the supercharges for some particular models such as the trigonometric models, 
and the Fateev-Zamolodchikov spin chain, and spin chain related to the supersymmetry t-J 
model, and a spin chain with mod — 3 conservation of the particle number. Moreover, we made 
a number of conjectures on the existence of zero-energy states and their properties. The most 
remarkable observations are certainly the sum rules for the singlets of the Fateev-Zamolodchikov 
chain, and a relation of its components to the weighted enumeration of alternating sign matrices. 

The present work leads to many open questions. We have seen that the solution of ([6| 
was instrumental in order to find spin chains with supersymmetry. Already at spin one, we 
observed the existence of several inequivalent solutions to the condition on the supercharges to 
be nilpotent. It is natural to ask if there is a way to classify all possible solutions for higher spin. 
Some of these solutions yield Hamiltonian densities for integrable/Bethe-ansatz-solvable spin 
chains. We observe this phenomenon in cases where additional constraints such as spin-reversal 
symmetry are imposed. It seems natural to ask if it is possible to set up a criterion for a 
supercharge to generate an integrable spin chain. Certainly, the analysis Reshetikhin's criterion 
for the two-site Hamiltonians would be an interesting starting point, yet more refined tools to 
prove integrability should be developed. 

For spin chains with lattice supersymmetry the possible existence of exact zero-energy states 
is certainly the most interesting phenomenon. Even though the conjugation argument provides 
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a tool to prove their presence in some cases, it does not cover all of them. For example, the 
XXZ spin chain with twisted boundary conditions cannot be analysed in this way as we do not 
know an ofF-critical extension which preserves the supersymmetry and the twist. Therefore, 
alternative existence proofs could put conjectures for twisted spin chains on solid grounds. The 
existence of supersymmetry singlets is intimately related to the cohomology ker QAr/imQ^v-i 
[55] . The number of zero-energy states coincides with its dimension. Looking at the list of 
known cohomology theories, the present study shows striking similarities with the structures 
used in the theory of cyclic cohomology [HT] . It would be interesting to establish a consistent 
connection with this theory. This might provide an algebraic tool to prove the existence of the 
zero-energy states. 

A symmetry is typically useful in order to determine quantities of physical interests such as 
ground-state correlation functions. Hence, it would be interesting to understand to what extend 
the supersymmetry presented here may simplify to the explicit evaluation of ground-state 
correlators. An example was provided in |25j where so-called scale- free expansions around 
trivially solvable points were pointed out. These allow to guess rather simple closed formulas for 
lattice correlation functions (such as the spin two-point correlations of the XYZ chain along the 
supersymmetric line at various distances |88p. For a direct application of the supersymmetry, a 
natural place to start would be infinitely long spin chains where the distinction between N and 
+ 1 sites is not an issue. Besides, the analysis of an infinite system N oo might prove to 
be interesting by itself: for example is it possible and/or meaningful to find a central extension 
of the lattice Af = (2,2) supersymmetry algebra? 

Finally, there is an obvious generalisation of the type of supersymmetry discussed in this 
article. We presented a theory based on a local operation q : V ^ V (>^V. It is natural to 
generalise it to q : V^" — ygim ^j^j^ arbitrary n < m. In this case, the supercharges add 
m — n sites to the system. The corresponding Hamiltonians will lead to interactions which 
involve in general more than just two neighbouring spins. 

Let us now turn to the specific models presented in this paper. We studied a special set of 
lattice supercharges, inspired from supersymmetric fermion models with exclusion rules. The 
sectors where the supersymmetry holds qualify as lattice-pendants of the conformal field theory 
Ramond sectors. Quite natural is the question: can we identify features of the Neveu-Schwarz 
sectors on the lattice? A glance at the spectrum of the spin— 1/2 and spin— 1 chains with 
periodic boundary conditions for small system sizes shows that they possess a single negative 
eigenvalue in the translationally-invariant sector if the number of sites is even. Hence this sector 
is a good candidate, and clearly deserves further investigation. Moreover, in this context the 
study of spectral flow (which will necessarily break the supersymmetry on the lattice), which 
may serve as a connection between the Ramond and Neveu-Schwarz sectors, appears natural 
[5U] . A further challenge is to construct in a systematic way the elliptic extensions of the 
trigonometric supercharges presented in this work, so that the lattice Af = (2, 2) supersymmetry 
algebra still holds. This was done in section |4] for the case of spin one. It seems natural to 
expect that these generate spin chains related to the spin— £/2 versions of the eight-vertex 
model at the special values rj = 'k/{£ + 2) of the crossing parameter, obtained through the fusion 
procedure. Indeed, the existence of a lattice supersymmetry at these points can be proved 
from the fusion equations, but it remains to work out the explicit supercharges. There are two 
possible routes to achieve this: either one uses the wave functions obtained from the algebraic 
Bethe ansatz, or one tries to solve the condition ([6| for a local supercharge q which changes 
the magnetisation by — (^-f 2)/2 mod 2. The off-critical models are intimately related to the 
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corresponding deformation in the Mi models. It was shown in [25 that for the spin— 1/2 the 
off-criticality on the spin chain side corresponds to stagger the fermion models with period 3 
on systems of length L — 3n. It is natural to extend this for higher spin to a staggering with 
period ^ + 2 on L = + 2)n sites. 

For the twisted Fateev-Zamolodchikov spin chain, the origin of the combinatorial features 
should be elucidated. A first step will be a proof of the sum rules, which seem to be a rather 
common features of both spin chains and fermion models with lattice supersymmetry |28j . 
Moreover, our findings hint at a simple eigenvalue for the transfer matrix of the 19- vertex 
model and its elliptic generalisation (a 21-vertex model [S5]) for any value of the crossing 
parameter and the elliptic nome. The exact diagonalisation of the transfer matrix for small 
system sizes confirms this conjecture. This allows to apply the machinery of the quantum 
Knizhnik-Zamolodchikov equation without any restriction to a special point in parameter space. 
This problem will be addressed in a future publication |90j . Eventually, as already mentioned 
in the introduction and further supported by the examples given in this work, there appears to 
be a relation between combinatorics and supersymmetry. This is not the first time that this 
phenomenon is observed as shows the example of the Veneziano-Wosiek model . It seems 
that this is not only a superficial coincidence and remains to be understood. 
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A Trigonometric supercharges 



In this appendix we explain a constructive way to arrive at (17 1. Our strategy is the following. 
We specialise first the matrix elements hi^r.s),{m,n) to the diagonal case (r, s) — (m,n), and 
construct a particular solution which is invariant under spin reversal. In the main text, we 
show that this solution satisfies indeed the spin-reversal symmetry even for off-diagonal matrix 
elements. We start thus by considering the matrix element bm,n = ^(Tra,n),(m,n)i given by 

h - n'^ ' ^ 

'Jm,n — "m+n+l,m 

\fc=0 fc=0 

First, set to = n = 0. We impose 6oo — and thus have 

2 _ \ " 2 




2 _ \ " 2 



fe=0 

Next, consider the equation 6o.f-„ — 6£.„. Using the previous result we find that 

£—n—l n— 1 

2 , t 2 _ t 2 , t 2 

k=0 k=0 
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Replacing now n ^ l~n and subtracting both equations, we find a useful symmetry property: 

^m,0 + "£+2-m,0 — ^1,0' m = 2, . . A. (35) 

Now impose 5^,1 — hi-rn,t.-i- For m = 0, . . . , — 2 it follows that 



^m+2,1 — 9 I ^ "f-m^A; T 2-^^l~l,k 





In particular, the case to = leads to a| = q = X]fc=o kl'^i thus to 

_^ / H—m — l m— 1 \ 

^m+2,1 = 2 I 5-/ ^i-m,k ~ ^ ^m.fe 1 + ^^2,0 ~ 2*^^'°' 
V A,-0 fc=0 / 

Again, we notice that the expression in brackets is antisymmetric under m — > ^ — to. This 
allows to derive the equation 

2 I 2 rj 2 2 



The left-hand side can be reduced with the help of the recursion relation (151, and thus 
everything can be rewritten in terms of q. After a little algebra, we find the non-linear 
recursion relation 

2 2 _i_/'2 2 V2 2n 2 fry 2 2 \ 

^m+2,0^m+l,0 + 1*^1,0 '^m+1,0 A'^l.O '^mfi) — '^1, 01^*^2,0 '^l,o)- 

We will now show that it admits a particular solution related to a linear recursion relation. To 
this end, set q — Zm/ Zm-i and normalise zq — 1. We find that 

Zm-l{Zm+2 ~ ^l^m+l + Z„i) + Zrn{Zrn+l ~ ^l^m + (2(^1 — Z2) — l)z„i_i) = 0. 

This relation holds certainly if we choose for Zm as solution of the equation 

Zm+2~ ZiZ,n+i + Zm^Q, with Zo=l,Zo = 0. 

This is the recursion relation for the Chebyshev polynomials of the second kind. We prefer to 
use g-integers and thus write 

= [to + 1], [n] = — j-, with q = e*''. 



The phase 77 is not arbitrary as we have to respect the symmetry relation ( 35 1 . After some 
algebra, one finds that this relation implies sin(£ -f 2)77 = 0. The unique choice for 77 to keep all 
the g positive is 77 = tt/{1 + 2). Hence, q takes the root-of-unity value 

Thus, we have [to -|- £ + 2] = —[to] = [—to]. Hence, we find 



The general coefficients am,k are obtained from the recursion relation (15 1. This completes the 
derivation of the constants ( 17 1 



42 



B Special polynomials 



Bazhanov and Mangazeev introduced in |5UI I51j a series of polynomials s„(2;), n ^1,, which 



appear as components in the ground-state vectors of the spin— 1/2 XYZ Hamiltonian (12 1, 
and sum rules (see also the related work of Razumov and Stroganov [61j ) . A sum rule for the 
square norm was recently proved by Zinn-Justin 52J. 

These polynomials appear to be not only off-critical generalisations of the combinatorial 
numbers found in the ground states of various variants of the spin— 1/2 XXZ chain at A = —1/2, 
but display also an interesting relation to classical integrability, namely to certain tau-function 
hierarchies associated to the Painleve VI equation. They are solutions to the non-linear 
difference-differential equation 

2z{z - l)(9z - if (In s„(z))" -|- 2(3z - lf{9z - l)(lns„(z))' 

+ 8(2n + i)2 -^"+i(^K-i(^) _ (4(3„ + i)(3„ + 2) + (9z - l)n{5n + 3)) = (36) 
s„(z)^ 

with so{z) = siiz) = 1. This equation can be understood as a result of special Backlund 
transformations for Painleve VI ^921 ES] • It is far from being evident that its solutions are 
polynomials for the given initial conditions. Various properties of the solutions to this equation 
were conjectured in f51|. For our purposes, the factorisation properties are relevant. In fact, 
the authors claim that for any fc G Z one has 

S2k+i{w'^) = S2k+i{0)Pk{w)pk{-w), pfc(O) = 1, (37) 

where the Pk{w) are polynomials of degree degpk{w) = fc(fc + l) with positive integer coefficients. 
We have the transformation symmetry 

6w\ I 1 ~ I 



Pk{w) = — pk 



The pk{wYs are the polynomials appearing in the zero-energy states of the twisted elliptic 
Fateev-Zamolodchikov chain at a; = ±1/2, y = ±C/2. 



C Counting momentum states 

In this appendix, we use the classical theory of necklace enumeration in order to count certain 
momentum states for a spin— i'/2 chain with N sites [M] . 

We consider thus a chain/necklace with N beads or equivalently a periodic lattice with 
N sites. Each site of this graph, labeled by some integer j = 1, . . . , A^, can be coloured 
with a colour nij = 0,1,...,^. Thus, a colouring/configuration is given by an A^-tuple 
/i = (mi, . . . ^niN-iTmN)- We denote by X the set of all possible configurations. There 
is a natural action of the cyclic group, generated by the cyclic shift T, on X through T • 
(toi, . . . , mAr_i, TTiAr) = {mpf,mi, . . . ,r7ijv_i). We will study the fixed points of this group 
action. This is motivated by the fact that the elements of X are labels of the Hilbert space of 
the spin chain. 
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Momentum states. To each /i = (mi, . . . , totv-i, mjv) G ^ we associate a basis vector 
l/i) — |toi, . . . , mAr_i, totv) in V^^ where V ~ C''^^ as in section [2] The translation operator 
on V'^^ acts according to T\^) = \T ■ ^) (we use the same notation for it acting on the vector 
space and the set X). Out of any basis state we build a momentum state 

N-l 

3=0 

where k = 2Trm/N, m = 0,1, . . . , N — 1 is the momentum. |^/;^) may be zero what implies that 
the configuration fi is incompatible with the given momentum k. In fact, for each configuration 
fi £ X there is a minimal positive integer r such that ■ fj, — fi. We call r the symmetry 
factor of /i. For the momentum state to be non-zero we need 

N 

r e with p 



'pd{N, to) 

Motivated by the main text, our aim is to count all momentum states with fc = for N odd, 
and k — TT for N even, with given total particle number 

N 

M = ^TO,. 

Chains of odd length. We start with N odd. A momentum state may be generated by 
two different configurations /i, /i' if they can be mapped onto each other through a suitable 
translation. This induces an equivalence relation on the set of configurations, and thus our 
task consists of counting the number vmn of equivalence classes (or equivalently the number of 
orbits in the set of configurations generated by the action of T) . 

We denote by X]^ C X the set of all configurations with given AI which are stable under 
the action of T™. According to Burnside's lemma, we have 

"MN = ^ E I^mI- (38) 

m=0 

Therefore, we need to evaluate |A"]J}|. The case to = is the simplest: we just need to count all 
possible configurations for given M. Let us thus consider to = 1: in this case all sites need to 
have equal colour because a single translation (cyclic shift) has only one cycle. Hence, Xl^ 
contains a single configuration of colour j if and only if M = jN for some integer j = 0,1, . . . ,£, 
in other words if N\M as long as M < iN. We abbreviate a^iM) = For to > 1 there 

various distinct cycles. It is a classical result that their number is given by gcd(A^, to), and 
therefore they have length r,„ = N/ gcd(iV, to). For each cycle, we have however the same 
problem as for to = 1. Hence, we find the general result 



\X 



Ml 



Mi,...,Mged(jv,m)=0 y 



/gcd(Ar,m) 

(39) 
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Wc notice a typical convolution structure. It is therefore useful to introduce the generating 
function An{z) for the numbers aniM) where z is the conjugate variable of M. We obtain 



An{z) = E ^n{M)i 



M 



Af=0 



1 - Z 



N 



We use this in ( 39 1 and find the generating series 



M=0 



1 — z^/ gcd{N,m) 



Performing the sum with respect to m leads to the generating function Jn{z) given in the main 
text for odd N . 



Chains of even length. Next, we treat the case of even N and momentum k ~ tt. Thus, we 
want to count the number of all representatives with even symmetry factor r. It seems natural 



to modify ( 38 1 in such a way that only fixed points of T™ with even m are retained. Thus, we 



attempt to discard all contributions with odd m and write ^/NJ2m=o^ I^m"I- "^^t, in this 
way we count some unwanted configurations. Indeed for m > 0, if we decompose 2m = 2^"^bm 
with hm an odd integer, and > 1 some integer, then the fixed points of T^™ contain the 
fixed points of T^™. These need to be subtracted. For to = we need to subtract the fixed 
points of T^"/2 . Thus we arrive at 




'^MN = ^ I > , I^M I - > , I for even N. 



The sum can be rearranged in a convenient way. To this end observe in (39 1 that \X^j\ depends 



on TO only through gcd(iV, to). In order to rearrange the sum we need the following lemma: 

Lemma C.l Let N be even and consider Cm — gcd(iV, for m — 1,2,..., N /2. We have 
Cm+N/2 — Cm- Furthermore, let p — then it follows that c„i+{p-i)/2 = gcd{N,2m — 1). 

Proof: The proof is entirely based on two simple identities. Let m and n be positive integers. 
(i) \i m > n then we may write gcd(m, n) = gcd(m — n, n); (ii) if n is odd then gcd(TO, n) = 
gcd(2TO,n). From the second identity it follows in particular that gcd(b„i,n) = gcd(TO,n) for 
odd n. Furthermore we find gcd(6m, n) = gcd(TO, 6„) = gcd(6m, 6„). 
For the first part, use 6jv = &Ar/2 and write 

c„i+N/2 = gcd(A^, b,n+N/2) = gcd(67v/2, ™ + N/2) = gcd(&jv/2, m) = gcd(6Ar, to) 
= gcd{N,bm) = Cm- 

For the second part, use 6jv — ^Af/2 = P and write 

Cm+(p-i)/2 = gcd (b^^p^,N^ = gcd ^TO + ^-y-;^JV^ = gcd(2TO + p - l,p) 
= gcd(2TO — l,p) = gcd(2TO — 1, bjsr) = gcd(2TO — 1, N)- 
This proves the claim. □ 
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Combining the implicit dependence of \X^\ c 
properties given by this lemma, we obtain 

N/2 N/2 
m=l m=l 

Putting everything together we thus find 

(N/2~l N/2 \ 

E i^II'vEi^mI 
m=0 m=l / 

Hence, we proved the extra factor (—1)™ for 
obtain the result given in the main text. 



Cm = gcd(A'', bm) with the periodicity and shift 



= i^E(-iri^Ml for even TV. 

m=0 

3ven A''. Computing the generating function we 
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